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ABSTRACT 


With the advent of inductively coupled plasma (ICP) 
emission spectroscopy, vast amounts of analytical information 
can be obtained simultaneously. However, only a small 
fraction of this information is actually utilized by the 
conventional detection systems such as direct reading type 
Spectrometers” when coupled’to the LCPs*in our “Laboratory, 
measurement systems based on photodiode array and Fourier 
transform spectrometers which are capable of simultaneously 
acquiring several hundred Angstroms of continuous spectral 
information from the ICP have been developed. Given such 
spectrochemical measurement systems, one then needs 
effective software and hardware systems to fully utilize 
the large amounts of spectral data available from these 
systems. Correlation based methods are found to be 
particularly-simple and-effective. 

In this thesis, theoretical and practical aspects of 
correlation techniques as applied to spectral data 
processing are discussed. Simply speaking, correlation 
analysis involves the use of a function (correlation mask) 
to extract desirable analytical information from raw 
experimental data. A unique correlation data point can be 
obtained by performing a point-to-point multiplication 
between the correlation mask and the raw spectral data at 


zero phase shift. The summation of these products represents 
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the consolidation of all the desirable spectral information 
into a Single point. With the proper choice of correlation 
masks, the relative and absolute magnitudes of this 
correlation data point indicate the qualitative and 
quantitative information for the analyte-of-interest. 

The correlation procedure is applied to analytical 
data obtained with the ICP-photodiode array spectrometer 
and the ICP-Fourier transform spectrometer systems. It 
is shown that this technique can be utilized to process 
complex spectral patterns for specific information via a 
software approach. This includes the identification 
and quantitation of sought-for analyte. Because of the 
Simplicity of this method, the possibility of implementing 
correlation analysis for automatic computer interpretation 
of spectral data is evaluated. 

Finally, the power of correlation techniques is 
then further illustrated with the development of a real 
time data processing system for interferometric signals. 

In this case, hardware implementation of correlation 
analysis is achieved with the aid of a low cost microcomputer 


and a high speed electronic multiplier-accumulator chip. 
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Vanadium analytical curves determined 

using nickel stripped vanadium binary 
cross-correlation masks. Thresholds (a) 20%, 
(b) 40%, (c) 60% and (da) 80%. 

Block diagram of the ICP-Fourier transform 
spectrometer system. 

Schematic diagram of the Michelson interfero- 
meter. 

Block diagram of the Michelson interferometer 
setup. 
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portional noise originated from the ICP. 
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a conventional argon ICP. 
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Analytical working curve for lithium obtained 
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Backeround of the argon ICP in ime visiple 
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Analytical working curve of calcium obtained 
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CHAPTER I 
Introduction 


A. An Ideal Spectrochemical Measurement System 

A typical spectrochemical measurement can be subdivided 
into three basic operations: an encoding step, a decoding 
step, and a processing step. At the encoding step, 
analytical information about the composition of a sample 
is encoded in the form of electromagnetic radiation. This 
part of the system is usually classified as the source and 
a common example is the use of flames in atomic emission, 
absorption or fluorescence spectrometry. At the decoding step, 
electromagnetic radiation carrying the analytical information 
is measured. This task is carried out by a spectrometer. 
In many cases, the processing step simply means the human 
interpretation of the spectral data obtained. With the 
advent of advanced electronic technology, a computer is 
usually incorporated into the system with the appropriate 
supporting data reduction software. It is then possible 
to obtain the desirable analytical results within a 
relatively short period of time. In addition, the computer 
may also-act) jas,a controller for the whole system as 
demanded by the user. 

Analytical spectroscopists have been in search of an 
“jdeal" system ever since the need for spectrochemical 


measurement arose. Although it is unlikely that an "ideal" 
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system will ever appear, it is worthwhile to look at some 
of the desirable characteristics of the "ideal" system to 
act as a standard towards which all "real" systems should 
strive. 

An ideal atomic spectrochemical measurement system 
Should have the capability of obtaining analytical 
information for all elements at all concentration levels 
and the data obtained should be a linear function of the 
concentration with no interference effects. In addition, 
accurate and precise results are desirable. The measurement 
system should also be applicable to all types of samples 
without much sample preparation. The instrument used to 
carry out the analysis should be compact, reliable, low 
cost and easy to maintain. With the ever increasing load 
of samples to be analysed, emphasis must also be placed on 
the necessity of the system to perform multielement 
analysis simultaneously and for final 
emaiveical results to. be. obtained in a short period, of 
time. This implies that efficient data processing methods 
must be employed. 

Among all the techniques that have been used to carry 
out the encoding step necessary for spectrochemical 
measurements, the observation of emitted radiation from 


free atoms generated from the sample offers the most 
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promise of providing an "ideal" system. Simultaneous multi- 
element analysis is possible with optical emission spectro- 
scopy in a relatively simple manner. The only requirement 
is an energetic source to generate excited free atoms 
and/or ions from a sample. In the next section, an almost 
ideal excitation source for atomic emission spectrometry 
(AES), the inductively coupled plasma (ICP), will be 


introduced and described. 


B. Atomic Emission Spectroscopy 

The principal excitation sources for AES are listed 
below: 

(a) Flames. 

(D) sparks. 

Vel bireee current (d.c,) arcs. 

(d) Plasmas. 

(e) Glow discharge lamps. 

(f) Hollow cathode discharges. 

(g) Lasers. 

Excitation sources based on glow discharge lamps, 
hollow cathode discharges and lasers have had somewhat 
limited usage in AES while d.c. arcs and sparks, together 
with the classical thermal excitation provided by high 


temperature flames such as nitrous oxide-acetylene, are 
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widely used in many applications. However, during the past 
decade, exceptional growth has occurred in the area of 
plasma emission spectroscopy. These energetic sources, 
especially the ICP are likely to be used for years to 

come as they approximate, in many ways, the ideal 
spectroscopic source. 

By definition, plasmas are gases in which a 
portion of the atoms or molecules is ionized. The plasma 
Sources that have been utilized in AES can be classified 
into: 

Cada Dec. tiphasmarge ts! 

(b) Capactively coupled microwave plasmas (CMP). 

(c) Microwave induced plasmas (MIP). 

(d) Inductively coupled plasma (ICP). 

An excellent description of these various kinds of 
plasma excitation source can be found in the article by 
Butler et al. (1). Comparison between these plasma sources 
along with a general review of plasma spectroscopy have 
also been presented by Greenfield et al. (2,3) and 
Boumans et al. (4). These references should be consulted 
for more details. 

The analytical performance of the ICP has been 
compared with the other plasma sources and with various 


established spectroscopic methods such as flame AES, flame 
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and furnace atomic absorption spectrometry, arc AES and 
x-ray fluorescence spectrometry by Boumans (5). The 
overall assessment was based on the analytical capability 
Lopperdorm, direct.solad enalysis; dinectwiliquids analysis; 
multielement analysis, trace analysis, precise analysis 
and. accurate. analysis.. He concluded that the ICP is: the 
method of choice although it suffers from the incapability 
of performing direct quantitative solid sample analysis. 
This clearly explains why the ICP occupies such an 
important role in AES at the present time. This can be 
verified by the degree of commercialization and ever 
increasing analytical applications of the ICP. 

Inductively Coupled Plasma: The ICP was developed from a 
low pressure discharge described by Bobat (6) in 1947 to 
the present atmospheric pressure gas discharge. The first 
ICP discharge sustained at atmospheric pressure was 
generated by Reed (7,8) in 1961. Reed used argon as the 
plasma gas and his system was applied to the manufacture 
of crystals from refractory materials. The easy 
vaporization of powders (A1,0.) injected into the discharge 
led to the suggestion that the ICP might have potential 
as an excitation source for AES. In 1964, the analytical 
application of the ICP discharge was first reported by 


Greenfield et al. (9). Shortly after, Fassel and Wendt (10) 
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published their own independent work on the ICP. Since then, 
an explosive growth in the use of ICP has occurred, 

A schematic diagram of the ICP discharge (11) is shown 
in Figure 1. It consists of an ICP torch which is composed 
of three ‘concentric-quartz tubes. The-torch is placed 
within a water cooled copper coil to which a radiofrequency 
(7-50 MHz, a typical value is 27.12 MHz) current is applied. 
The power output of the radiofrequency generator is usually 
in the range of one to three kilowatts, though some ICP 
systems operate at power levels as high as ten kilowatts. 
me ICP is initiated witngtneraid tof a Tesla coil. It 
provides a seed of electrons which interact with the 
oscillating magnetic field (see Figure 1) produced by the 
high frequency current to further ionize the plasma support 
gas (normally argon) thereby forming a stable discharge. 

The temperature of the plasma discharge in the central 

core 1s about 10,000 °K, Three flows of gas, usually argon, 
are introduced to the torch. The plasma gas (or coolant gas) 
is introduced tangentially at about 15 liters per minute (lpm) 
and is intended to sustain the discharge in an annular 

shape and to prevent the outer quartz tube from melting. 

The auxiliary gas (0-1 lpm) is used to raise or lower the 
plasma discharge and in some cases, this flow of gas is 


not necessary. Samples, in the form of fine aerosol are 
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Schematic diagram of an ICP discharge. 
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introduced into the discharge with the aerosol carrier or 

Mebulizer gas through ‘the central tube of the ICP torch: 

The optimal emission region of the plasma is normally 

around 15 mm above the copper coil depending however on 

the operating conditions and the sought-for analyte (12). 
some of the desirable characteristics offered by 

the ICP as an encoding device in a spectrochemical 

measurement system are: good sensitivity and low detection 

limits for most of the elements, high precision and accuracy 

in the results obtained, linear response with respect to 

analyte concentration of up to five orders of magnitude, 

and low interference effects. In a sense, the ICP closely 

approximates the "ideal" spectroscopic source other 

than the fact that at the present time it only produces 

Sood. analytical results for liquid samples and it is 

relatively expensive. Also as a result of the high excitation 

@nerey avalilablevin the ICP, spectral interferences can be 

a serious problem. However, the ICP does provide the 

spectroscopist with an excellent excitation source for 

spectrochemical measurements. The next problem is how to 

extract and utilize the vast amount of analytical information 

available from the ICP. In other words, what kind of 

spectrometer should be used to decode the analytical 


spectrochemical information generated by the ICP? 
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C. Multielement Analysis Spectrometers 

ineorder Vo iudily utilize the spectral information 
available from the ICP, wide continuous spectral coverage 
by the spectrometer is highly desirable. Furthermore, 
Simultaneous detection capability is necessary for efficient 
spectral analysis. This is also essential for background 
correction which is a major consideration when using an 
ICP. A true background or spectral correction procedure 
can only be implemented if simultaneous spectral coverage 
Can be erie eae 

In general, spectrometers commonly used in AES can be 
Subdivided into two main types: 

(a) Dispersive type. 

(b) Non-dispersive type. 

Dispersive type spectrometers have been 
conveniently grouped into four categories by Boumans (13). 
The list is reproduced below because it gives a general 
idea ane various kinds of dispersive spectrometer that 
are available. 

(a) Polychromators, one-dimensional dispersion with 

fixed optical components: 
(i) large polychromators (direct readers) with 
multiple exit slit and multiple detector 


(photomultiplier tube, PMI) or multiple 


oad waiter! iadrogte. vo ala ant ‘oa 
sap revao [Setorgs suaun ines lw pa ait oe ' 
atom 191 Teed i § alone baa al ‘te ape 2 
sig Rte 1dt yuesEees bah sd pe tenemes insaiee soe J seu re 
Bighichan sed sinikd itrieaeha 3 oaie ag bia vatayeana , 
' . ata z haisa * ai soi a 
axibeoats roles one ae a. ru mtoad 9 
lion banat SN Be cvosnetlumie . 3 Kategertae os 


af dso STA wi fess ie Seuninm's ne. cae 
| ae Ba etna ih ai 
ee sen ey ui 


. (ef) semaie® ay +95 
anaes 2BVEN: $2. Ses 
al bagi’ oth sige kb 


% 


fit tw Roteiayeth teem artsinib = 


ie ‘ 


sotaeteh Seazatas = Fite % 


(b) 


Gc) 


(d) 


10 


exit slit and single detector (PMI) or a 
Single movable exit slit and a single detector 
(PMT) ; 

(ii) large or small polychromators with exit slits 
replaced by multiple arrays of photodiodes, 
each consisting of a few photodiodes for 
measurements of the line peak and the back- 
ground. 

Echelle spectrometers, two dimensional dispersion 

with fixed optical components: 

(wes) Multiple exat slit and single detector (PMT) 
design, “spinning encoding disk"; 

(11) multiple exit slit and multiple detector 
design, cassette with fixed apertures, mirrors 
and PMTs; 

(iii) with television camera tube (and image 
converter; 

(iv) with image dissector. 

Spectrographs, one dimensional dispersion with 

photographic Breet on completed with: 

( i) manual microphotometer; 

(ii) recording microphotometer; 

(iii) computerized automatic microphotometer. 


Monochromators, one-dimensional dispersion with 
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mechanically movable optical component(s) for 
wavelength adjustment: 
(i) single exit slit and single detector (PMT) 
design with wavelength control as follows: 
-manual 
-scanning 
-slewing and programmable with precision 
mechanical device 
-slewing and programmable with computer or 
microprocessor; 
Gix) without. exit) slatinwithetelevision camerador 
self-scanning linear photodiode array. 

Among the dispersive type spectrometers, there is no 
doubt.about the-.capability of spectrographs in covering °a 
wide spectral region, simultaneously and continuously. In 
fact, spectrographs, which use photographic detection can 
record thousands of spectral lines within a short period 
of time. However, the major drawback in the use of 
spectrographs is that time consuming readout procedures 
are required. Although computerized automatic microphoto- 
meters such as those developed by Witner et al. (14) and 
by Walthall (15) have appeared to minimize the problem, 
good quantitative results are still difficult to obtain. 


These disadvantages can be overcome by the use of detectors 
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such as photomultiplier tubes (PMTs). With the use of these 
detectors, electronic readout can be instantaneous. Linear 
response over 5-6 orders of magnitude with good quantitative 
results can be obtained with PMTs. However, with the use 

of PMIs which are essentially single channel detectors, wide 
spectral coverage capability is lost. In order to combine 
the advantages offered by PMI's and photographic detectors, 
direct reading type spectrometers were introduced. In this 
case, up to sixty channels of PMTs can be configured in 

the focal plane of a polychromator for simultaneous detection 
of spectral information. The major disadvantages offered 
by the use of this type of spectrometer are: non-trivial 
alignment procedures for each channel, limited PMT channels 
because Of phys. cal-Limrcavcons In “therfocalsplane or the 
polychromator, the instrument is bulky and expensive, 
adert 1s Gltirrcurt to change vo new channels = These, 
combined with the fact that continuous spectral coverage 

is not possible for background correction are some of the 
reasons for the appearance of slew scan type spectrometers. 
For slew scan type spectrometers, a single channel detector 
is used to scan through the entire spectral region. This 

is achieved either by moving the optical components of 

a monochromator or by moving the detector along the 


focal curve cf a polychromator. In this way, continuous 
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spectral coverage can be achieved. However, precise 
alignment of the detector at a particular spectral wave- 
length for analysis is not easy even with the help 

of a computer. Furthermore, the capability of simultaneous 
detection of signals is lost. Much time is spent in scann- 
ing through the entire spectral region in search of the 
sought-for analytical information. 

So, the problem of combining simultaneous, continuous 
and wide spectral detection with good quantitative measure- 
ment characteristics into a single spectrometer has puzzled 
spectroscopists. for, along time... With. the arrival of 
electronic image detectors such as self-scanning photodiode 
arrays, image dissectors and television camera type detectors, 
an extensive investigation on the feasibility of these 
kinds of detectors for analytical applications has begun. 
These detectors are essentially the electronic versions 
of photographic detectors. They were developed to combine 
the advantages offered by both PMTs and photographic detectors, 
although: the goaljis. not. quite reached, at, the present, time. 
These devices are capable of simultaneous coverage of a 
limited spectral region with instantaneous electronic 
readout. Within the last few years, developments in this 
area have produced image detectors such as the intensified 


photodiode array with a sensitivity comparable to PMTs. 
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The major drawbacks of the image detectors are high cost 
and limited spectral coverage. Despite this, these electronic 
devices have found numerous applications in the development 
Gia  taeal = spectrometer. Of particular anterest is the 
replacement of single channel detectors such as PMTs with 
amage detectors. They can be used in the direct reading 
type spectrometers such as those developed by Boumans 
(16,17) or slew scan type spectrometers developed by 
Peclack ey al, (lo, 9) Fry (20), Yates et:al. (21): and 
Furuta (22). With these configurations, simultaneous 
coverage within a narrow spectral region for background 
Correction is possible: In order to extend the spectral 
BOvVora26 Ol these cCetecvors, ine use of Bcheile spectro- 
meters which produce two dimensional spectra combined with 
Ewe use of two dimensional image detectors has been 
Reporved. Lhic, contieuration was) first. reported by 
Margoshes (23) in 1970 and since then, several versions 
have been developed (24-27). With the tremendous progress 
in electronic technology, it is likely that the potential 
of the image sensors in the development of an "ideal" 
spectrometer is not fully developed. 

All the spectrometers mentioned above belong to the 
dispersive catagory. Another approach is to dispense 


with the dispersive system completely and make use of the 
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multiplex technique. With this technique, spectral informa- 
tion is encoded so that a single detector can be used to 
simultaneously measure a wide spectral region. The most 
common examples are Fourier transform spectroscopy (FTS) 

and Hadamard transform spectroscopy (HTS). In the case of 
HTS; a dispersive element is still used, In fact it is very 
Similar to a direct reading spectrometer and is listed under 
the dispersive type spectrometers by Boumans. The only 
difference with the HTS is that a single detector (PMT) is 
used and discrete spectral information from various exit 
Slits of a polychromator are multiplexed into the single 
detector. Since a’ multiplex advantage is not realized in 

the UV-visible-near IR region where detector noise is not 
theslimiting» factors; HTS does not. offer a better.alternative 
to the conventional dispersive type spectrometers when 
coupled to the ICP. However, in the case of FITS, wide 
continuous spectral information can be obtained simultaneously. 
The use of a FT spectrometer to carry out atomic emission 
spectrochemical measurement has been discussed in detail by 
Horlick, Hall and Yuen (28). It appears that the major 
drawback of the FITS when used in conjunction with the ICP is 
the dynamic range problem. This limitation is due to the 
multiplex property of the technique. Every single data 


point contains information covering the whole spectral 
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region and must stay within the scale of the measurement 
electronics system in order to obtain proper analytical 
information concerning the sample. The dynamic range 
problem exists if the sample gives a very intense emission 
at a particular wavelength. This intense emission will 
Memb tewne amplitication of the signal and emission from 
weak spectral lines is difficult to detect. With the use of 
mre ICY, the dynamic range problem due to strong background 
emission can essentially render some spectral regions 
Boalivyiically useless as reported by Hall (33). .This can 
mow be resolved with the use of a mixed gas ICP. Further 
Srecucslon ot thre problem can be found,.in. Chapter VL. 

in conchusion,. there 1s)mo definite answer as to which 
design is the "ideal" spectrometer. However, one thing that 
Gan be said fe that simultaneous and continuous spectral 
coverage is a major feature of an "ideal" instrument. 

In our laboratory, spectrometers based on photodiode 
arrays and Fourier transform techniques capable of 
Simultaneous multielement analysis have been developed. 

The photodiode array spectrometer can simultaneously 
measure 1024 spectral resolution elements and the Fourier 
transform spectrometer, 4096 spectral resolution elements. 
Compared to multichannel polychromators which may have 


20 to 50 channels the measurement capability of these new 
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systems represents a major increase in capability. Given 
this new level of spectral information measurement 
capability one then needs effective software and hardware 
systems to fully utilize such large amounts of spectral 
data. A basic step in the processing problem is the 
G@evectron oO: the Spectral information Characteristic of the 
chemical species. In atomic emission spectrometry this 
typically involves the detection and quantitation of a 
specific combination of spectral peaks in a complex 
multiline spectrum. 

In this thesis, correlation techiques as applied to 
the data processing step of the spectrochemical measurement 
systems based on ICP-photodiode array spectrometer and 
ICP-Fourier transform spectrometer will be presented. 
Correlation techniques were chosen because they are simple 


and highly effective in processing complex data arrays. 


D. Correlation Based Signal Processing 

Correlation techniques have long been used to measure 
and process signals in physics, chemistry and engineering. 
An early work concerning the application of correlation 
analysis to the detection of periodic communication signals 
was reported by Lee et al. (63). This topic was reviewed 


in some detail in 1960 (64). Correlation techniques, again 
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applied to communications, have also been described by 


Lange (65). Many of the applications of correlation techniques 


to chemical data have been reviewed and presented by 
Horlick and Hietfje (38). Interested readers should refer 
to this review and the references cited for full details. 

Correlation techniques can be effectively applied to 
spectrochemical data. In particular, the Signal detection 
aspects of correlation analysis look promising for achieving 
computerized automatic detection of spectral information. 
The power of correlation techniques for signal detection is 
best illustrated with an example. 

Three atomic emission spectra are shown in Figure 2 
for zinc, cadmium and boron all measured in the same 
spectral region (~220 nm) with the photodiode array spect- 
rometer described in Chapter III. The wavelength axis has 
been left off on purpose to emphasize the fact that 
knowledge of it is unnecessary for spectral identification 
using cross-correlation techniques. Let us choose zinc as 
the sought-for spectral information. The complete cross- 
correlation patterns for zinc-zinc, zinc-cadmium, zinc-boron 
are shown in Figure 3. The exact computation of these 
cross-correlation functions will be described in the 


next chapter. 


Cross-correlation patterns can be intuitively evaluated 
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FIGURE 2. (a) Zinc, (bd) cadmium and (c) boron spectra. 
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FIGURE 3. Cross-correlation patterns for the (a) zinc, (bd) 


cadmium and (c) boron spectra with the zinc spectrum. 
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by envisioning the correlating pattern, the zinc spectrum 
in this case, to be slowly translated across the other 
spectral patterns. The value of the cross-correlation 
function at any particular displacement is the sum of the 
products of the two functions at that displacement, in a 
sense, their mutual areas. A unique point is the point 

of zero displacement, i.e. when the wavelength axes of the 
two spectral patterns exactly coincide. This point is often 
designated the t = 0 point. A large relative value of the 
cross-correlation ane os at t = 0 indicates a high 

degree of similarity between the two patterns. This is 
Clearly indicated in Figure 3 by the cross-correlation 
pattern of zinc with itself. Note the relatively large 
central value and the symmetry characteristic of an auto- 
correlation pattern (Figure 3a). Thus the relative magnitude 
of the t = 0 point can be used to indicate detection of a 
spectral pattern. Note that the other cross-correlation 
patterns have no distinct maximum at t = O although the zinc- 
cadmium cross-correlation pattern does contain a peak near 

t = O indicating that two lines of zinc and cadmium must 
exist at essentially the same wavelength. Such is the case as 
zinc has a line at 213 nm and cadmium a line at 214 nm in 
this region. 


The absolute magnitude of the t = 0 point is also 
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important in that it represents the consolidation into a 
Single value of all the JERR Me information about the 
sought-for constituent. Thus rather than quantifying zinc 
on the basis of a single spectral intensity, the t = 0 value 
represents the summation of the quantitative information 
m@rom all spectral lines, thus more completely utilizing 

the available spectral information. 

Correlation techniques can be further extended to 
process complex patterns for specific analytical information. 
This task includes the identification and quantitation of 
sought-for analyte with the proper use of cross-correlation 
functions. Further details and examples of correlation 
techniques as applied to atomic emission spectrometry will 
Peepresented ever sin Chapver lil. First,” the theory 
of correlation analysis will be reviewed. In addition 
the mathematical computation of the cross-correlation 
functions will be presented. A software approach of 
implementing cross-correlation techniques for processing 
spectral data obtained from an ICP-photodiode 
array spectrometer system and an ICP-Fourier transform 
spectrometer system will be described in Chapter III and IV 
respectively. A real time data processing system for 
interferometric signals utilizing cross-correlation 
techniques implemented in a hardware approach Wale be 


described in Chapter V. 
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CHAPTER II 


Theoretical and Practical 


Aspects of Correlation Analysis 


The basis of correlation techniques, as applied to 
chemistry, has been well documented in the literature (34-39). 
In this chapter, correlation theory based on these references 
will be summarized and presented. For a detailed description 
Or pie Theory and its applications to chemical data 
measurement, the review by Horlick and Hieftje (38) should 
be consulted. As will be shown later, the full correlation 
Operation can be carried out with the aid of Fourier 
woansiormavion. Ihe practical computation, of the correlation 
function via the Fourier transformation route will be 
illustrated with particular attention paid to the manipula- 


Lon. i0f .imput, and joutpuw idata, arrays. 


A. What is Correlation? 

Simply stated, correlation analysis provides information 
about the coherence within a signal or between two signals. 
The correlation function of two signals is obtained by evalu- 
ating the time averaged or integrated product of the two 
signals as a function of their relative displacement. 
Mathematically, correlation can be expressed as: 

nate 


Cap(t) = lim = a(t) b(t+t) dt (a8) 
Tro 2T J -T 
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where Cap(®) is the correlation function between the two 
Signals a(t) and b(t) and t is their relative displacement. 
Pee Signals acanebe a function.of essentially any variable, 
for example, wavelength, retardation, frequency, accelerating 
Voilvare, time te. sthus,,.if a-and Db .are considered to .be 
PubcCoLOUSsOL ime selLne .correlation.function Cob will be 
related sto .,and,plotted, against the relative time delay 
between, the two.signals. 

In most situations, correlation is implemented on 
degitized.signals.iFor-digitized signals, .the calculation 
of the correlation function can be expressed by the 


following summation: 
Captnahiic= D a(t)b(ttnat) H=OFAG2 Ope. MELA GHR) 


The signals can only be displaced some integral number of 
ppeasampling interval, Av. Thus, the displacement’ nwt is 
ecoulvaleny. vO 1t In Equation i. 

There are two basic correlation operations: auto- 
correlation and cross-correlation. These, together with a 
special kind of correlation operation, convolution, will 
be treated individually later. First, the relationship 
between correlation and Fourier transformation must be 


defined. 
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B. Correlation and Fourier Transformation 

The ready availability of digital computers has 
facilitated software approaches to correlation. In 
particular, the fast Fourier transform (FFI) algorithm 
(40) has essentially revolutionized the extent to which 
correlation techniques can be implemented on computers 
and indeed many correlation based techniques are more 
commonly referred to as*Fourier transform techniques. 
The-power and range of applicability of some of these 
techniques are generally well known and documented in 
several books (41-45). In addition to providing a 
convenient route for the implementation of correlation, 
a knowledge of Fourier transforms also aids the utilization 
and understanding of correlation operations. An important 
theorem concerning correlation states that correlation of 
two waveforms is equivalent to multiplication of their 
Fourier transforms followed by inverse Fourier transformation 
Gn tne produel.* sSchematrcalty, this Sequencercan be 


represented in the following way: 


a(t) te b(t) = Cab(Tt) C152) 
FOURIER INVERSE 
FOURIER 
TRANSFORMATION TRANSFORMATION 
= Can(f) 
NEw, x B(f) AB( (ay) 
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The asterisk in Equation iii is merely a common 
Shorthand way to represent correlation; Equation iii is 
therefore-..dentical ‘to “Equation i. ‘Notice that ‘in 
Equations 111i and iv, cross-correlation is shown to 
involve a multiplication process in the Fourier domain. 
simple das this aspect of correlation ‘might seem, it is 
of paramount importance. A large fraction of software 
and hardware methods for obtaining the correlation function 
rely on the Fourier transform method because of its speed 
and relative simplicity. 

several versions of the FFT abound in the literature 
and in program libraries and utilization of these sub- 
routines may, at first, seem trivial. However, subtleties 
and pit falls exist with respect to data pretreatment, 
data post-treatment, inverse Fourier transformation and 
manipulation of real and imaginary arrays that can cause 
considerable grief to the uninitiated. Thus, the utilization 
and recycling of FFT algorithms will be discussed before 
correlation analysis is illustrated. 

Computer Software: All calculations were carried out 
using a FORTRAN FFT program which is based on a Cooley- 
Tukey FFT algorithm (40) first suggested by Gentleman and 
Sande (43). The complete program listing with flowchart 


is included in Appendix III. This program can be operated 
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with or without a sorting routine which will be explained 
later. The maximum number of points that can be processed 
with the DEC PDP 11/10 minicomputer running under an RT-11 
operating system is 4096 due to the in core memory size 
(32 K) limitation. Detailed documentation can be found 
along with the program. All the figures are 1024 points 
unless otherwise specified. They were displayed and 
plotted with a VT-11 screen and a Zeta incremental 
plotter. 

Uiulize ton Of -rastarourver Transtorm Algorgtnms: AVL 
FFT algorithms require real and imaginary input data arrays. 
However, in the vast majority of situations, the pros- 
Bective user has a simple array of data to process; for 
example, aldigimtized optical spectrum. Tnus the user is 
immediately faced with the problem of how to fill up the 
real and imaginary"imput larrayseoftthe FFT with the array 
of data. With most FFT algorithms, it is recommended 
taswetneudatayarvay poe placed in the real) inputearray and 
Phoabathen|macinaryyimput array be filled with zeros. This 
is illustrated in the following scheme.with R and I 


representing the real and imaginary arrays. 


R 
FFT 
(1) 
DATA 


I (ZEROS) ai 
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A data array which is a simple line spectrum generated by 
wien eOmpurer Vs used for illustration and is shown in 
Figure 4. If this data array is placed in the real input 
arrayj;ot the FRI and zeros,in-the imaginary input array 
Unesresulting real and imaginary outputs, after a 1024- 
pol; Frit, fappear as shown in Figure 5. 

An alzernate method of filling {the real and imaginary 
input arrays with data to be transformed is to put the 
CMCMmDOlINUGh 2.6.1, 1, Jou, 7s.» Of the data array into the 
fen snoutiarray andethe even points, i.é.; 2, 4, 6, Bs... 
OB @ne data array into the imaginary input array. This 


scheme is outlined below: 


R (odd points) FFT R 
DATA | > (2) 
WITH SORT 


I (even points) 


When this method is used the FFI must be followed by a 
Sorming routine in, order to generate|the comrect veal and 
imaeinary outputs (46) . 

This method provides for very efficient use of the 
FFT algorithm as N data points can be transformed using 
an N/2 FFT. However, this approach is normally carried 
out so as to keep the transform size identical to the data 


array. This is easily accomplished by filling out each 
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FIGURE 4. A simple line spectrum used as the input data 


array for the FFT implementation. 
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input data array to N with zeros after placing the N/2 
odd and even points in the real and imaginary inputs. 

ihe results of this approach are illustrated in 
Pigure 6. The same input data as that illustrated in 
Figure 4 was used. The 512 odd points were placed in 
the real input array and the 512 even points in the 
imaginary array. Both input arrays were then filled out 
to 1024 with zeros and a 1024-point FFT calculated 
followed by the sort routine. 

Both Schemes 1 and 2 are equally valid approaches to 
setting up the real and imaginary input arrays of the 
FFT algorithm starting with a simple data array. However,- 
scheme 2 issusually preferred because the real and 
imaginary OuopuLS are intuitively easier to anterpret 
in that the end of each output array represents the 
Nyduist frequency rather than the center of the array 
as for Scheme 1. For example, the FFI outputs resulting 
from Scheme 2 implementation for a single peak spectrum 
eiownein Pigure 7 are illustrated in Figure co. (The 
real and imaginary FFI outputs are damped cosine and 
Sine waves. The functional form and extent of the damping 
are directly related to the position of the peak along 
the horizontal axis (wavelength axis in the case of optical 


spectrum) relative to the origin. The amplitude output 
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FIGURE 7. A single peak spectrum. 
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FIGURE 8. (a) Real, (b) imaginary and (c) amplitude output 
arrays for the single peak spectrum with Scheme 


2 FFT implementation. 
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array Pere + 12)/2; is shown in Figure 8(c). It 
represents the amplitude of the Fourier components making 
up the original data, the single peak spectrum. In all 
cases, these outputs, using Scheme 1, would be complicated 
by information at each end of the output arrays. This is 
illustrated in Figure 9. These are the real and imaginary 
output arrays of the single peak spectrum obtained 

using Scheme 1 implementation. 

Inverse Fourier Transformation: With a great many 
Fourier transform data processing methods, data is 
transformed, manipulated in the transform domain, and then 
re-transformed (inverse Fourier transformation) back to 
she original “domainv@ Ihe basrseformthese methods involves 
phesitact that convolutions, cross-correlations, and 
auto-correlations can all be performed as multiplications 
in the Fourier domain as pointed out earlier in this 
chapter. 

A number of subtleties exist in implementing inverse 
Fourier transformations that may give many a novice 
considerable frustration. Several alternate routes are 
possible depending both on the original method of Fourier 
transformation (Scheme 1 or 2) and on the nature of the 
processing carried out in the Fourier domain. The main 


confusion centers around the exact manipulation of the 


os) 


ft .{o)8 onmget ai awoda el ‘el ot oF | 

siden etascoqnes Sekwot eat Yo sbettique edt eam 

tia «at -eeitoods aa0q sfania ried ated Lai 

beteo hinges sd bilvow yk amade? pitew Pt ucreene wa 
ei 22? oye Fingitre gt? 3s fae tege wa wed 


rivabanal bis fas ote sta spe? .g etuslt ae 


: 
soniesdo om ehoads daeq aiantis eds ts 2¥ 
_ febetatoonstgnt £4 
Yam teotg 2 TY same egpee tS eos. gat 2 - 
gi sted ,abodoen iitnedpirtn. aio avotenere § 
nade bie jlamoh wpetanesy sit al heraliy Lome 8 om P: 


f 


ey Hoge (acivamnotagatt setusd% eexovel) boone ) 2 
seviovel sbhbotten Awsntt <0? otasd ont cb aioe {er 
. bres , a0 ibale thn 3015 dota kovema pad? tom 
eteltecviigition as ‘bomb rueq 24 Siege arel tater 

ain? mf 1sifiss Fug hak iog oe BisKeh sonia 
sareval giitmeualqe? mi. fads ee tS “rodeniet A 

epivon « yam svtq, gan Feat encipsorrotenett xt U9 

ats detu0r odentsits Laxeved .nottererr? afsaxehbede 
wizyol to bodren Laehgixe eaF a6 Mtod yrtbasgsh aldzesog 
eres | 


Mage ieddele 


Git to eutan ert? no Rite 4S co £ seatot) motdnenotenet 
aien sai? -oteno® abit ait a¥ tivo boteses galeassey 7 
ent to noltalivginet dane eas Brvons ree Cote 


36 


“UOTIe{USWSTAUT Tq JT ewWayog 


YZIM UNITZoeds Yyeed aTSuts 9uy aoz skessze qndyno ArTeutseut (1) pue Tesey (e) *6 


(Q) (2) 


08 t6— ee a 


2 ec— 
(@1X) SCAT TdWy 


Q@ BE 


( 


QB 26 


G8 BS' 


4ynolt 


00 OST—= 


08 86— 


@ BE 00 @E— 
NOX) 30NLT TdWy 


( 


88 06 


0 @c' 


ot 


+h 


. pybrremer eae 


TRI MI PP 


- 
cal 


ayuey6 Leap ober 


ie 


tot 


— 


n 


hig hie 


OAT 


abe 

a 
<8 

' 


TRSs 


*. 


2) Tee 


= 


7 


a. 


——w <a he wh “a een 


~t3094 - 


eo . Cin 
‘ i‘ 7 a = Bs Fi 
REPT it. Weer) 


36 


explicit real and imaginary arrays that are generated 
after the first PFT. 

The standard route for recycling the FFT when 
eeheme 1 16 used jiofyset up the imput arrays for the first 


transform is shown below: 
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ine real output array from the first FFI as used as the 
real snpur Tor athe second FFT and the negative of the 
imaginary output array from the first FFT as usedas the 
imaginary input for the second FFI. The real, imaginary 
and amplitude output arrays of the second FFI using Scheme 


3 are shown in Figure 10. The original data is the line 
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FIGURE 10. (a) Real, (b) imaginary and (c) amplitude 
output arrays for Scheme 3 inverse FFT 


implementation. 
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Sie) 


spectrum illustrated in Figure 4. In this case the real 
and the amplitude output arrays both contain the original 
data. 

When Scheme 2 (even-odd point method) is used to 
set up the original input arrays, recycling is carried 


out as outlined below: 
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zeros to the N value of the transform. The real, imaginary 
and amplitude output arrays of the second FFT are shown in 
Figure 11, again for the case in which the original data 

is the line spectrum shown in Figure 4. Now only the 

real output array contains the original data. 

In addition to Scheme 4, several other recycling 
implementations when Scheme 2 is used as the primary FFT 
method have been tried. These results will be presented 
in Appendix I as they may aldvothers an Sorcingsowe Pr. 
problems and options. 

Though the same pattern of the original spectrum 
can be obtained using different versions of inverse FFT, 
one has to watch out for the change in amplitude. In 
general, the amplitude of the resulting spectrum is 
increased by a factor which is equal to half of the number 
of relevant data pointsinputted for the inverse FFT. For 
example, in Scheme 3 the amplitude of the final spectrum 
is'a factor of 1024 higher than the original spectrum 
since there are 2 x 1024 (real and imaginary) data points 
inputted for the inverse FFT. In Scheme 4, the amplitude 
MerOuueDvec taeu0r Of Jie because Onity tone VO2h "data array 
has been used. Thus, amplitude adjustment is necessary 
whenever an inverse FFT is implemented. This applies to 


correlation operations if the Fourier transformation 
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route is chosen. 


C. Cross-correlation 

A cross-correlation function is obtained if the two 
signals, a(t) and b(t) in Equation i (see Section A) 
are different. This*operation-is used to check for 
Similarities between the two signals. 

Full cross-correlation in the Fourier domain requires 
complete systematic utilization of the real and imaginary 
arrays of the FFT algorithms. The basic general scheme 
that is used for full cross-correlation operations is shown 
in Figure 12. The two functions to be correlated (function 
one and function two) are processed with the sorting FFT 
algorithm. Each will produce real and imaginary arrays; 
RL, 11 and R2;"i2) "ay comphex mathematical multiplication 
is then performed. The product now consists of real, 
CaaiR2 + 1112) “and imaginary, (R211.- R112) (parts. These 
are inputted into the appropriate real and imaginary input 
arrays of the non-sorting FFT (algorithm. iInethe case of 
the imaginary array, the negative is taken. The 
resultant correlation function will then be available in 
the real output array resulting from the non-sorting FFT 
algorithm. The whole operation is illustrated in Figure 


13 with two simple functions. The functions shown in 
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(a) Function one and (b) function two for cross- 
correlation operation. (c) Cross-correlation 
function as calculated by the procedure outlined 


iirigure 12. 
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Figures 13(a) and 13(b) are simple spectral patterns. 

Cross-correlation patterns can be intuitively 
evaluated by envisioning one of the functions, for example 
function two to be slowly translated across function one. 
pee value of the, cross-correlation function at any 
particular displacement is the sum of the products of 
the two functions at that displacement, in a sense, their 
mutual area. In this particular example, when the two 
functions are cross-correlated, the resulting correlation 
function obtained in the real output array as shown in 
Figure 13(c) is out of sequence when plotting against a 
scale running from a negative displacement (-t) toa 
positive displacement (+t). The actual correlation pattern 
of the two functions is shown in Figure 14. It is obtained 
after proper rearrangement of the data array. Rearrangement 
simply means the interchange of the upper 512 points with 
the lower 512 points of the function which contains the 
correlation information shown in Figure 13(c). 

At this point, the distinction between a(t) * b(t) and 
b(t) * a(t) correlation operations should be made. This 
is best illustrated by the above mentioned example. If 
the two functions are interchanged in the correlation 
operation, the resultant correlation function will then 


be the mirror image of the one obtained previously. This 
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is shown in Figure 15. Thus, care must be exercised 
when performing full cross-correlation operations as to 


which signal among the two is displaced. 


D. Auto-correlation 

Auto-correélation is cross-correlation of a function 
wath itseit. If the two signals, a(t) and -pét)-an Equation 
eeare identical, the funetion obtained is’ called mn autoe 
correlation function. It is used to check whether coherence 
exists within a signal. It can be represented by the 


following equation: 


Uae Pah 
Caa(t) = lim — a(t)a(t+t)dat (v) 
T +00 al =p 


maere Caa(t) is the auto-correlation function fori the 
Signal a(t) and t is the relative displacement. 


In the case of digitized signals, the equation becomes: 


Caa(nat) = }) a(t)a(tenat) T0512 aeons ee) 
t 


Auto-correlation can be implemented by the scheme 
shown in Figure 16. It is somewhat simpler than cross- 
Perrelation in that only two FFT's need to be carmied out 


and there is no imaginary cross product. The multiplication 
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product which contains only real terms will generate the 
auto-correlation function in the real output array 
Pesulting irom the non—sorting FIT algorithm. “in this 
case, the imaginary input array for the non-sorting FFT 
algorithm is filled with zeros.- The rearrangement of the 
eorrelation Tunczion output array is necessary as in the 
case of cross-correlation. An even simpler route is 
possible for auto-correlation using only the sorting FFT 
elgorithm and thig 2s illustrated in Figure 17.4 A simple 
spectral pattern and its corresponding auto-correlation 
function are shown in Figure 18. The full auto-correlation 
memctlon as depiered in Figure 19 can be obtained simply 
by reflecting the resulting function in Figure 18(b) about 


the origin. 


E. Convolution 

A number of the aspects of correlation that have been 
introduced are perhaps more familiar to the reader under 
the term convolut ore However, convolution, an operation 
which occurs during the generation or measurement of all 
signals, can be considered to be merely a special kind of 
correlation process. Mathematically, convolution can be 


expressed ass 
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If Equations vii and i are compared, their only 
eeirerence, is observed tobe a. minus sign ian front: ofthe 
for, the. Dotunction~ The errect oflthisiminus sign” is to 
reverse the b(t) function (from left to right) on the time 
Exiles betore at is multiplied by a(t). Thetrest of the shirt— 
ome, Multiplying and averaging procedures are identical to 
those employed in correlation. Therefore, correlation and 
convolution are identical, except that in convolution, one 
ef the signals is first t-reversed. 

The reason for this can be best understood with the 
aid of an example.:- Consider the effect. of anvinstrument 
response function on a signal. A peak input signal anda 
low pass filter response function typical of amplifiers 
and recorders are shown in Figure 20. By convention Fi 
tradition both the signal and the response function are 
expressed with the time increasing to the right. However, 


if we imagine the input signal entering the instrument it 
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is apparent that the left side of the signal would enter 
first. Thus, in order to use the correlation operation to 
assess the effect of the instrument on the signal either 
the signal or the response function must be reversed from 
Pett to right before evaluation. It is conventional to 
Mevcvee Unie response function. Thus, convolution is cross- 
correlation with reversal of the response function before 
evaluating the correlation function. The convolution 
function of the peak signal and the low pass filter response 
function obtained using the same flow chart (Figure 12) as 
in the case of cross=correlation is illustrated in Figure 
zi. The response function is reversed before the full cross- 
correlation operation. 

From these considerations, convolution can be seenito 
be merely a special kind of correlation and can often be 
implemented and utilized in similar ways. In fact, in 
those many instances in which the response function is 
symmetrical, its inversion produces no change, so that 


correlation and convolution provide identical results. 


men information av vt = 0 Point of the Cross-corréelation 
Function 
In the previous chapter, the power of correlation 
analysis for spectral identification of relatively simple 


spectra has been illustrated. The t = 0 points of the 
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cross-correlation functions contain quantitative and 
qualitative analytical information corresponding to the 
Sought-for analyte. In this section, another example is 
given to further illustrate that in many cases, full cross- 
correlation functions need not be evaluated. As long as the 
two functions that undergo cross-correlation operation are 
Eoepnase, the cross-correlation function at the’ t = 0 

point contains all the relevant analytical information 
required. This is best illustrated with a more practical 
erampie: the qualitative and quantitative analysis ofa 
multielement sample spectrum. A spectrum containing 

Caqgmium and Zine is shown %in Figure 22a,This, together with 
the zinc, cadmium and boron spectra shown in Figure 2 were 
measured with the ICP-photodiode array spectrometer system 
in the same spectral region (~220 nm). The experimental 
details in obtaining these spectra will be given in the 
mext chapter. 

The spectrum containing cadmium and zinc information 
is cross-correlated with the pure Zinc, cadmium and boron 
spectra individually. At this point, the term "correlation 
mask" should be introduced. It simply represents the function 
used to mask out desirable analytical information from a 
sample spectrum by correlation techniques. In this particular 


example, cross-correlation masks are essentially the zinc, 
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FIGURE 22. Sample spectrum containing cadmium and zinc. 
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cadmium and boron spectra which are used to extract zinc, 
cadmium and boron information respectively. This terminology 
Will be frequently used throughout this thesis. The results 
of the cross-correlation operations is shown in Figure 23. 
The relatively high magnitude of the t = 0 points correspond- 
ing to zinc and cadmium correlation masks clearly indicates 
mhe presence of Zinc and cadmiuméspectral information inthe 
jo Plecar ne tie caceeot. Doron, noidistinct nigh sc = 0 point 

aS observed indicating the unlikelihood that boron is present 
in the sample. The steel nee magnitude of the t = 0 point 
will provide the quantitative information corresponding 

to the specific correlation mask used. As long as the 
amplitude of the cross-correlation mask remains constant, 

the magnitude of the t = 0 point will be linearly 
proportional to the quantity of the sought-for information 
present-in the signal. 

Ghissexample 1s given here merely 10 stress the 
capability of spectral identification and quantitation 
using the cross-correlation technique. The t = 0 point can 
be utilized as opposed to the full cross-correlation 
operation as long as the correlation mask and the signal 
are in phase. Full cross-correlation or scanning between 
the signal and the mask is not necessary. In the subsequent 


chapters, correlation based data processing utilizing the 
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FIGURE 23. Correlation patterns for the sample spectrum with 


the (a) zinc, (b) cadmium and (c) boron spectra. 
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t = O value of the full cross-correlation function as applied 
to multielement atomic emission signals originating from 


ean ICP source will be illustrated and discussed. 
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CHAPTER lit 


Correlation Based Data Processing for 


the ICP-Photodiode Array Spectrometer System 


A. Introduction 

Cross-correlation techniques have been effectively 
applied to the computerized automatic detection of spectral 
intormation (36,48,49). The basic approach as outlined in 
the last chapter involves cross-correlation of the raw 
spectral signal with a noise free mask of the sought-for 
Srectrakbipattiermiad Irnsthivs ychapter themcoricept sof Ycross- 
correlation methods for spectral pattern detection as 
applied to atomic emission spectra measured with a 1024- 
element photodrode array spectrometer coupled to an ICP 
will be outlined and then several analog and binary 
spectral masks will be evaluated as to their effectiveness 
Hor spectralorsienalydetection. yin particular special 
etLenta ongisi paid tho theiveffectiveness of cross-correlation 
techniques in cases where spectral interference is severe. 
This is important because when electronic image detectors 
(50-55) are coupled to spectrometers, spectral resolution 
is often sacrificed in order to achieve reasonable 
wavelength coverage. 

The cross-correlation techniques will be illustrated 
for the determination of nickel or/and vanadium. This 


system was chosen because the determination of both elements 
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is of prime importance for the petroleum industry (56-58). 
Besides, both nickel and vanadium yield complex spectra 
which result in several spectral overlaps. An evaluation 
Oppineset tectiveness,of cormrelation,.analysis.for the 
automatic processing of spectrochemical data obtained 


by an image sensor can thus be made. 


B. The ICP-Photodiode Array Spectrometer System 

A commercially available radiofrequency (27.12 MHz) 
inductively coupled plasma which is capable of delivering 
a maximum power of 5 KW was used as the source for all the 
atomic emission spectra described in this chapter. The 
plasma was imaged onto the entrance slit of the photodiode 
SGhaywspocrromet en asic. 4UV gerade -quantz.~lens of. .cen 
centimeters focal length. The whole experimental setup 
Wash 10Gated on nanuoppLeal.rai) ped similargin concept.to 
the one suggested by Walters (59). 

The 1024-element self scanning photodiode array 
spectrometer has been described in the literature (18,75). 
A brief description of the spectrometer is included below 
for completeness. For full details, references cited 
should be consulted. 

A 1024-element photodiode array was placed horizontally 


along the exit focal curve of a monochromator. The 
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monochromator is of the Czerny-Turner design and has a 
focal length of 35 centimeters. The standard diffraction 
grating provides approximately 500 A of spectral coverage 
Hor the array detector. The physical length of the arnay 
is approximately one inch while the individual photodiode 
measures 0.00) inch wide by 0.0147 inch gall. During data 
PeGuisition, the gpnotcdiode array was cooled to about: 220°C 
Minn the use of a Peltuer cooler. “<This significantly 
reduces dark current of the array and enables the 
anvegratyon time*to be: varied fromimilliseconds: Go Seconds 
depending on the magnitude of the signal. The photodiode 
array spectrometer was coupled to a PDP 11/10 minicomputer. 
The complete system with the supporting software has been 
recently described by Hull (60) and will not be repeated 
henge. 

The analog signal obtained from the spectrometer was 
electronically filtered and amplified with a differential 
amplifier prior to analog-to-digital conversion. The 
analog-to-digital converter was contained in the PDP 11 
Laboratory Peripheral System. The equipment used for this 
study is summarized in Table l. 

The atomic emission spectra of nickel and vanadium 
were obtained at an ICP power level of 1.5 KW with the 
coolant gas flowing at about 15 liters per minute (lpm). 


The auxiliary gas flow was completely shut off. A 
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conventional cross-flow nebulizer was utilized with the 
gas flowing at about 1 lpm. The observation height of 
the plasma was approximately 15 mm above the load coil. 
The spectrometer was set so that it covered approximately 
from 275 to 325 nanometers. All the spectra obtained were 


averaged over ten seconds unless otherwise specified. 


C. Simple Cross-correlation Masks 
theitigst step in the generation of'a eross—correlation 
Mesk tor;the detection of specific spectral information is 
the measurement of a spectrum of the sought-for element. 
The atomic emission spectrum of vanadium is shown at the 
bottom Of Froure 24, *The waveTeneth axis waseleft off. 
When utilizing cross-correlation techniques it is not 
necessary to define the wavelength axis, the specificity 
is built into the measured vanadium cross-correlation mask. 
This complete spectrum can be used as a cross- 
Peer tation mask to detect vanadium in subsequent 
measurements. Only the t = O point need be evaluated 
which can be generated by multiplying this spectral 
intensity cross-correlation mask and subsequent spectra 
together point-by-point and adding up the products. Such 
a procedure can be used to establish an analytical curve 
plotting the magnitude of the t = O cross-correlation point 


versus concentration. Such an analytical curve is shown 
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in Figure 25. Alternate masks may be generated using 
pereshold devels. Several such levels are indicated in 
Figure 24 superimposed on the vanadium spectrum. Now, 

gn the generation of the t = 0 value; cross-correlation 
products.are only evaluated when the vanadium mask value 
exceeds the specified threshold. With the 80% threshold 
value the measurement reduces to essentially a conventional 
Single peak intensity measurement as might normally be 
Garriedgout in setting up an analytical curves The 
analytical curves for various threshold levels are also 
Shown in Figure 25. As less and less vanadium information 
Pomusecd=cor establish the analytical curve the sensitivity, 
as indicated by the slope of the analytical curve, is 
reduced. Thus’ the full cross-correlation mask provides 
the most sensitive measurement of vanadium. 

The full cross-correlation mask also provides the 
most precise measurement of vanadium. Ten vanadium 
spectra at the 10 ppm level were measured and the t = 0 
cross-correlation points evaluated using the masks with 
various thresholds. The percent relative standard 
Gevietion (ed) of the ten t = .0 values ranged ‘from 11.6% 
for the 80% threshold mark to 8.3% for the full (0% thres- 


hold) cross-correlation mask. 


An interesting alternative to the "spectral intensity" 


atiay bovaiseag sd) yam aitaaan oderrres iA “— te 

ei betso lent 35 adevek dosne farveves -otorae & 
wou rere cttheney aay: qe bewogais 3g a 
si soeaatin 90 =#g072 alas a =v ont to nol 
siiev Mon abéGne siisy E ss al fos eeLevs vite | ‘ ras 
ptoviesipft we ont aE eae naitinega # 
iraaines 6 I LS tease Sei a trsmartisodenn » 


ais +eyia: gs Let tiugtead ke ae nt iobe ae 
Oaks 9%e ebares b Ladieeaibeld euadesy fot — 

me tfemrate 2b i there a49f ( end ang. Ba 3 
erky Lr ketise sat) sygia a oe 
=p jay tie sctigtana. bas, Vio “qote any yi £ 
nak ive teens aa ae ber one wins * 
“pig SDALEM, +8) c is 7 
nis sebivotg ceES oa nei 
ou Lea eget + sg a: r 
Os Feit bigs beet: pot hee be on ak 
si 208 a¥asu of anise bate | 


ee 


‘ VOny, - ibe Mis 
byebnsse 7 7 


see 30) Ly? ‘sat Sam 


i 
SyTianstat eu ne “ia i, wnt ba 
tacetat Leatsuage" seit oy ae 


3.47 


3.68 


2.69 


2.30 


1.91 


CORRELATION DATA 


1.52 


100.00 128. 02 199.00 168.08 162.88 200.00 
PPM VANADIUM 


PIGURE 25. Vanadium analytical curves determined using 
spectral intensity cross-correlation masks. 
Threshold (a) 0%, (b) 20%, (c) 40%, (d) 60% 
and (e) 80%. 
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masks depicted at the bottom of Figure 24 are the binary 
masks shown at the top of Figure 24. A binary cross- 
correlation mask is generated by first setting a threshold 
and then every spectral intensity greater than or equal 

to the threshold value is set equal to binary 1 and every 
mete less than the threshold is set equal to binary O. 
Binary masks for thresholds of 20%, 40%, 60% and 80% are 
snown at the toprof Figure 24. 

Sevenalwadvantages accrue to the utilization of binary 
Cross-eorrelation masks. The evaluation of the vt = 0 point 
rs Moumeoreideranly simpler as the multiplication step is 
not required. The t ='0 cross-=correlation of the binary 
mask and a measured vanadium spectrum amounts to a gated 
gnvegratton Ote specific spectral intensities. Analytical 
curves established using the binary masks are shown in 
Figure 26 and are comparable to those shown in Figure 25. 
Also the precision of the t = 0 value is essentially the 
same for the binary and intensity masks. At a threshold 
Value of 20%. the percent rsd of the: tengp = O values for 
the 10 ppm vanadium spectra was 8.7% for the spectral 
intensity mask and 8.5% for the corresponding binary mask. 

In addition to allowing simpler evaluation of the 
cross-correlation values the binary masks occupy con- 


siderably less storage space in a computing system than 
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FIGURE 26. Vanadium analytical curves determined using 


binary cross-correlation masks. Thresholds 


(a) 20%, (b) 40%, (c) 60% and (d) 80%. 
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the spectral intensity masks. This could facilitate 
implementation of cross-correlation methods with simple 


ROM (Read Only Memory) based microprocessor systems. Also 


integrated circuit subsystems are now available specifically 


designed to carry out the cross-correlation between a 
stored binary pattern and an input analog signal (61). 
Utilization of these subsystems would allow essentially 
real time evaluation of t = O cross-correlation values. 

The analogous spectral intensity and binary cross- 
Correlation masks for;nickel in this exaclh.same spectral 
region are shown in Figure 27. They can be used for the 
determination of nickel in a manner analogous to that just 
discussed for vanadium. 

These cross-correlation masks are good for spectral 
meentitication and quantification for relatively, simple 
Snadysis.) When analyvical information is required for 
samples of complex nature, spectral interference is a 
major problem. In the next section cross-correlation masks 
with Duilt—-in specificity for a particular analytical 
problem, determination of vanadium or nickel in a mixture 


Of two, Willi be illustrated. 


D. Cross-correlation Masks for Complex Spectra 


Consider the problem of generating cross-correlation 
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masks for the determination of either nickel or vanadium 

in the presence of each other. The spectrum of a solution 
containing both nickel and vanadium is shown in Figure 28. 
several lines overlap, the most serious being the overlap 
of the most intense line of nickel in this region (310.1 nm) 
with the second most intense vanadium line (310.2 nm). 

A spectral intensity cross-correlation mask specific 
for vanadium in the presence of nickel can be constructed 
by stripping a vanadium spectrum with a nickel spectrum. 
This is easily accomplished by subtracting a nickel 
spectrum from a vanadium spectrum as shown in Figure 29. 
much an operation is very easy to Seer aR in a precise 
manner using the computer coupled photodiode array 
epeetrometer. Depending on the degree of spectral overlap, 
ume amount of Stripping ‘can be easily adjusted. In this 
example, the stripping was carried out using a nickel 
spectrum whose maximum intensity was twice that of the 
maximum intensity in the vanadium spectrum. This was 
accomplished by computer software on stored spectra. The 
remaining features in the vanadium spectrum (above a 
threshold of zero) can be used to set up specific spectral 
intensity and binary cross-correlation masks that will 
detect vanadium in the presence of nickel and will not 


respond to any nickel information. These masks are shown 
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in Figure 30. The analogous procedure can be used to set 
Up cross-Ccorrelation masks for the detection of nickelflin 
the presence of vanadium and the results are shown in 
Figures 31 and 32. Thus with the measurement capability 
of the photodiode array spectrometer it is relatively easy 
t@ Se€t up, in an essentially automatic fashion, cross- 
Gerrelation masks for the analysis of mixtures. | Again, 
note that knowledge of the wavelength axis or even where 
mines Overlap is notarequired. The stripping procedure 
automatically removes the overlapped spectral intensities. 
hep et = 0 value ‘of the’ cross-correlation function 
was evaluated for spectra obtained of solutions containing 
a constant amount of nickel (1000 ppm) and increasing 
amounts of vanadium (100 to 200 ppm) using cross- 
correlation masks generated from a nickel spectrum and 
a vanadium stripped nickel spectrum. The responses 
of these spectral intensity correlation masks 
are shown in Figure 33 and 34 respectively with 
various threshold levels. In all cases,, pure nickel 
spectral intensity cross-correlation masks do extract 
undesirable vanadium information from the signal as 
indicated by the slope of the graphs in Figure 33. The 
vanadium stripped nickel cross-correlation masks do not 


respond to vanadium information as illustrated in Figure 
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FIGURE 33. Detection of nickel in the presence of vanadium 
using nickel spectral intensity cross-correlation 
masks. Thresholds (a) 0%, (b) 20%, (c) 40%, 
(d) 60% and (e) 80%. 
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FIGURE. 3) .«.Detectaon.of nickel inthe presence of vanadium 
using vanadium stripped nickel spectral intensity 


cross-correlation masks. Thresholds (a) 0%, (b) 


20%,.(c).40%, (d).60%.and (e). 802. 
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34. Similar results are obtained for binary cross- 
correlation masks at different threshold levels, these 
eeesillustrated ad bigures 35 and 36. Linear working 
curves for these cross-correlation masks are illustrated 
for the determination of vanadium in the presence of 
nickel as shown in Figures 37 and 38. 

At this point, one can envision correlation masks 
designed to implement more complex decisions than simply 
the detection of a single element. This can be achieved 
by Simultaneously assessing several elemental species and 
uaeir relative -eoncentrations. In this way, correlation 
masks with specificity built-in can be generated for a 
specific analytical problem. 

One has to keep in mind that the correlation masks 
generated as mentioned above can overcome the problem 
of spectral interferences, they do not compensate for the 
true chemical interferences arising due to the matrix of 


the sample. 


E. Conclusions 

From the fact that the performance of binary cor- 
relation masks is comparable to the spectral intensity 
correlation masks in all aspects, a preliminary investigation 


to generate binary cross-correlation masks from a library 
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Detection of nickel in the presence of vanadium 


PUGURE 95, 


using nickel binary cross-correlation masks. 
Thresholds (a) 20%, (b) 40%, (c) 60% ana 


(ad) 80%. 
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FIGURE 36. Detection of nickel in the presence of vanadium 
using vanadium stripped nickel binary cross- 
correlation masks. Thresholds (a) 20%, (b) 40%, 
(c) 60% and (da) 80%. 
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FIGURE 37. Vanadium analytical curves determined using nickel 
stripped vanadium spectral intensity cross- 
correlation masks. Thresholds (a) 0%, (b) 20%, 

(c) 40%, (d) 60% and (e) 80%. 
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FIGURE 38. Vanadium analytical curves determined using 
nickel stripped vanadium binary cross- 
correlation masks. Thresholds (a) 20%, (b) 40%, 
(c) 60% and (da) 80%. 
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of spectral line wavelengths by computer software was 
attempted. This was accomplished by calibrating the 
spectral window covered by the photodiode array with 

some standards that emit known spectral lines within the 
particular spectral window. Appropriate correlation 

masks were then generated for data processing. This has 
not been very successful because of difficulties in exactly 
calibrating the photodiodes to corresponding wavelengths. 
This is due to the fact that the dispersion of the 
monochromator varies slightly with wavelength. If this 
would have been successful, it would partially solve the 
problem of limited spectral coverage by the array detector. 
The monochromator can be slew scanned and calibrated 
before cross-correlation masks for specific analytical 
problem will be generated for data processing. 

Based on the results obtained in this study, 
experimentally generated cross-correlation masks can be 
very powerful. However, more spectral coverage by the 
image detector is necessary in order to fully extend the 
advantage of correlation analysis. Image sensor based 
spectrometers utilizing area arrays and Echelle spectrometers 


are now being developed that could eliminate this limitation. 
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The fact that the knowledge of wavelength axis is not 
necessary may help spectroscopists to further extend the 
analytical capability of this system. Spectral information 
can be packed into the detector without regard to the 
wavelength axis to fully utilize the size limited array 
detector (62). The task of extracting the specific 
epectral.pattern, for,.each element); can; then be handled. by 
cross-correlation techniques even with severe spectral 
overlaps. Without cross-correlation analysis, 
human interpretation would be fairly difficult since all 
the spectral information does not necessary follow 
the wavelength axis. 

Gorrehationtanalysasedoes not-Jdamitsatselixuto,be 
used as an-effective data reduction tool for the photodiode 
array spectrometer system, it can also be implemented in 
ausimilarfeoncept forsether: spectrometers. basedson 
different image detectors. In the next two chapters, 
cross-correlation techniques are utilized to process 
spectral data arrays in the form of interferograms generated 
by Fourier transform spectrometers, With this kind of 
spectrometer, wide spectral coverage and precise wavelength 
calibration can be obtained in contrast to the dispersive- 
based detector system such as the photodiode array. Thus, 


the potential of correlation techniques can be further extended. 
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CHAPTER IV 


Correlation Based Data Processing 


for the I1CP-Fourier Transform Spectrometer System 


Pe inuroduction 

Simultaneous wavelength detection capability is one 
Ot the chief advantages offered by a Fourier transform 
spectromever when coupled to an ICP. A wast amount of 
analytical information is available within a relatively 
Short period of time. Thus, effective data processing 
is of prime importance when the system is utilized in 
multielement atomic emission spectroscopy. This 
Statement is also true in other fields of spectrochemical 
measurement utilizing an FT spectrometer as the detection 
system, such as in the area of infrared spectroscopy 
where the FI spectrometer finds most of its application. 

The time domain signal (i.e. the interferogram) 
obtained from the spectrometer is normally processed 
through the fast Fourier transform (FFI) mathematical 
algorithm described earlier in Chapter II. This step is 
necessary before human interpretation of the spectral 
data is possible. The FFT, when utilized has imposed 
certain limitations on the efficiency of the data 
processing routine, in addition to the requirement of a 
certain computer memory size for the processing of the 


spectral data. In the case of a large number of data 
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points for high resolution work, a large computer memory 
is desirable. 

Hardware development in computer technology, such 
as the array processor, has tremendously reduced FT 
computation time required by the computer. But, these 
devices are usually very expensive. In the case of 
long interferograms, techniques such as decimation in 
time (76-78) where interferograms are processed in 
segments rather than as a whole giant block of data can 
be employed. However, software development was found 
to be non-trivial (33). The laboratory computer can 
also be interfaced to a large central computer for efficient 
processing of large blocks of data points. In these cases, 
efficient data processing will be hindered by the time 
used to transfer data between devices. In view of the 
above comments, a simpler, alternate method to efficiently 
process the spectrochemical information obtained from the 
ICP-FT spectrometer system is sought. Capability of the 
data processing method to perform real time analysis is 
of prime consideration. The adaptation of the method to 
automatic interpretation of spectral data is also 
important. Correlation is the answer. 

Correlation techniques have been applied to process 


interferograms (33,79,80). In this study, the feasibility 
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of the cross-—correlation function at the te= 0 point as the 
data processor for multielement atomic emission inter- 
ferograms will be investigated. The advantage of the 
precise wavenumber axis provided by the FI spectrometer 
will be utilized. This enables cross-correlation masks 
to be generated by the software of a small laboratory 
Computer wor; extraction of analytiical information 
directly from the time domain signal without going 
through the FFT process. As mentioned above, the 
possibility of real time data processing and automatic 
interpretation of spectral data will be assessed. 
Ingorden to tully-appreciate this study, a general 
discussion on the use of the FT spectrometer as a 
detectionssystem for the] I1GP_is included?in @his chapter. 
Particular attention is paid to the dynamic range problem 


imposed by the use of the FI spectrometer. 


B. The Experimental System 

A block diagram of the experimental setup is 
illustrated in Figure 39. The inductively coupled plasma, 
Michelson interferometer and all the related optics 
were mounted on an optical rail bed similar to the one 
designed by Walters (59,81). The rail bed was in turn 


placed on a vibration isolation table which was air-floated 


28 


ait a8, Intog 0 <7 Sah Teme e san meee keer pee 
Ta Tc rotamer agers ; Ys Sitlem wt 30 


steoetort ae he 

ashton suodkig adage caae: th won 

att  svous eS HaESOONe, TR ost ty 
sitemodve bets satabaadaay st88. at tees ‘te Vin 
“baazozal goebay ave Seuropigh soit: sderte 
fstaisg. a a yeeTS td ote seins: abot oo f co fel 


" 
A 


Tet ner TOR ges i‘ keivelondh ait ae 
venget®. air ni Pababhe ag, i faa apraya ai 
TaALYosE ayQins 1 + SA Rectyp alt BF vttter oa 
OS Vee om v8 


ee 


fe 
. 
/ 


m. 


aw >, capgmi-bepig ge 39 einsdd took 44 
pemaadt belaies ne es Saat i sg 


$40. oAT 
ory fs 


batsoll-tis Saw qom 


94 


‘wa4sks Joqoulorzj0eds wWaoysuedy Jaetanog-got auy fo WetseTp MOOT *66.0N0DTA 


A\ddns 
wniyjeH 
/ue601NIN 


J S9jUu0I9e13 
=] Juewesnsee,; 


Ol/tt ddd 


030 
seyndwogd 


40}8J0U05) 
du 


JOX1IN 
/810\@W MO} 4 


$9/U02}90)3 
jOuyuod 


1s9ndO 


3088} ON-OH xog yew sue} sue} 10JOWOsEjE}UI 


| peg 118y 
| yuowuBbliy ddl BulBew; Bujyyeuyjjoo uOs|/®YUIIW 


-oui* 


Az2 


= 
a 


Puen GE, 


id ae 


» « he cf 
BS dy 
a 


= 


on 


* 
a. ay, 


‘ie 


whenever experimental data were being taken. This system 
provides precise alignment with the aid of a He-Ne 
laser and the excellent stability that is required for the 
operation of the interferometer. This setup is advantageous 
because alignment between the plasma and the interferometer, 
together with the internal optical alignment within 
the interferometer is very critical especially when it 
is being used in the ultraviolet region because of the 
Short wavelengths that are involved. 

The Michelson interferometer with its control 
electronics used in this study has been well documented 
an the leterature (26, 32,393). (For the latest modifications 
of the interferometer system, refer to Appendix II. A 
Schematic diagram ofthe interrerometer design and a bilock 
diagram of how the interferometer operates are included for 
completeness in Figure 40 and 41. The inductively 
coupled plasma is commercially available and is capable 
of a maximum power output of three kilowatts. The 
operating conditions for the ICP are similar to those 
described in Chapter III. In some cases, gases other 
than argon such as nitrogen or helium have been introduced 
into the coolant stream of the plasma. In these cases, 
analytical data were obtained under a mixed gas plasma 


condition. Mixed gas plasmas were initiated using the 
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FIGURE 41. Block diagram of the Michelson interferometer 


setup. 
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conventional argon and were then slowly crossed over 
to the mixed gas condition with the aid of a flowmeter. 
Further discussion on the use of the mixed gas plasma 
WillJappear later in this chapter. 

A small region of the plasma was spatially selected 
by a 10 cm focal length quartz imaging lens and a 5 mm 
Circular aperture. The radiation passed by the aperture 
ime ycOllimated by another quartz lens of 12 mm focal 
Penmezth beiore it reached the interferometer. This optical 
link between the ICP and the interferometer is not the only 
evernativen. Currently, an all mirror based optical 
interface is being considered. 

Different detectors can be used in conjunction with 
the interferometer, depending on the spectral region of 
maverest. A summary of the difterent detectors thar 
have been used to cover the near IR-visible-UV spectral 
Merion 16 ineluded in Table 2. Due to the fact that the 
present beamsplitter in the interferometer does not 
operate below 250 nm, the solar blind PMT detector is 
sometimes replaced by a combination of the visible PMT 
and a solar blind optical filter (Oriel 5786) which has 
a spectral response from 250 nm to 360 nm. This enables 
a gain in sensitivity of several orders of magnitude. 


However, with UV studies below 250 nm, the combination 
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TABLE 2. Detectors used in the ICP-Fourier transform 


spectrometer system. 


Detector Spectral Response Peak Response (nm) 
Solar Blind PMT MOO 320 220 


(Hamamatsu R166) 


UV-Vis PMT 200 - 650 300 
(RCA 1P28) 
Silicon Diode 400 - 1100 850 


(Electro-nuclear 
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is not feasible. 

The analog signal was fed into a current amplifier, 
electronic filter and a low noise amplifier before the 
analog to digital conversion was performed by the 
Laboratory Peripheral System that was interfaced to a 
PDP 11/10 minicomputer. All the interferograms were 
averaged over twenty five scans and output through an 
incremental plotter. Signals in the spectral- domain 
obtained by the FFT of the interferograms were obtained 
using the software developed by Hall (33). Detailed 
descriptions of the individual programs can be found in 
reference 33. Sources of the equipment that have been 
used in this system are summarized in Table 3. 

Two spectral systems were used for this correlation 
study. In one system the multielement signal obtained 
from a solution of nickel (1000 ppm), vanadium (100 ppm), 
cobalt (2500 ppm) and iron (1000 ppm) was studied. 

This signal was measured using the 1P28 PMT combined 

with the solar blind filter and a conventional argon 

ICP. The second system studied was the multielement 

spectrum of the alkali elements (1 ppm sodium, 2 ppm 

potassium, 0.1 ppm lithium, 10 ppm rubidium and 5 ppm 
cesium). This signal was measured using a silicon 


diode detector and a mixed gas nitrogen Ji arzon dor. 
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C. Use of a Fourier Transform Spectrometer as a Simultaneous 
Detection System for the ICP 
(a) General Considerations 

When an FT spectrometer which is essentially a 
Michelson interferometer is utilized as a measurement 
system for atomic emission spectroscopy, it offers various 
advantages over the conventional detection systems. Most 
of all, its wide and continuous spectral coverage is 
unmatched by any other spectrometer. In addition to 
this, high resolution can be achieved in a compact 
system, the resolution function can be controlled by the 
use of apodization techniques (32) and accurate and 
precise wavelength calibration is inherent in the nature 
of the instrumentation. All these, with the computerization 
of the system, provide a simplicity and flexibility that 
cannot be matched by other systems. 

Like most systems, the FT spectrometer is not 
without problems. One of the problems being aliasing, 
which refers to the phenomenon of undersampling of 
the signal. This problem exists especially when the 
spectrometer is being utilized in the ultraviolet-visible 
near IR region. It can certainly impose spurious 
information when dealing with broadband information. 


However, with atomic emission spectrochemical measurements, 
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line spectra are normally obtained. Yuen and Horlick 
(32,82) showed that alaising can be used advantageously if 
it is carefully considered and manipulated. Another 
problem deals with the degradation of the signal-to-noise 
ratio when the spectrometer is used in the UV-visible 
region. However, Signal-to-noise ratio is often not the 
overwhelming consideration when carrying out a spectro- 
chemical measurement. 

The necessity of carrying out the FFT data processing 
could be considered a drawback of the system because of 
the requirement of computer and software overhead. 

With the remarkable progress in computer technology, this 
problem will certainly be phased out. Above all, a 
Simpler, alternate data processing approach based on 
correlation techniques can be utilized. This correlation 
technique can be implemented in real time even in the case 
of long interferograms. Furthermore, it can be utilized 
in such a way that automatic computer interpretation 

of spectral data is possible. Full details on this 
technique will be described in this and the following 
chapter. 

By far, the most severe drawback that limits the 
practical use of the ICP-FT spectrometer system for real 


sample analysis is the presence of a dynamic range 
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problem. This limitation is due to the multiplex 
property of the spectrometer. Every single data point 
contains information covering the whole spectral response 
of the system and the signal must stay within the scale 
of the data acquisition system in order to obtain proper 
analytical information concerning the sample. The 
dynamic range problem is realized if the sample gives 
avery intense emission at a particular wavelength. This 
intense emission will thus limit the amplification of 

the signal and emission at some other weak spectral 

lines will thus be very difficult to detect. In a sense, 
the detection of weak signals is limited by the bit 
resolution of the analog-to-digital converter being 

used. This is why it is called the dynamic range 
limitation. This problem can be partially Overcome 

to aicertain extent by the use of optical .or electronic 
Gelters. ,,Howevenr, when the FI spectrometer 2s coupled 

to the ICP, a serious problem arises because of the 
presence of very intense plasma background (line 

and continum) extending from the UV to near IR region. 
Detection capability of the system may be severly limited 
by this background. Use of optical or electron cet ters 
to resolve the problem is out of question in certain 


spectral region because of the enormous background 
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and the close proximity of the spectral lines. Another 


important consideration is the "visibility" of the noise 


present in the ICP by the FT spectrometer. This is 


very important with respect to the use of the spectrometer 


as a simultaneous multielement detection system for the 
Powe Jittle bat of Gnsight into this 1s essential 
as the way the noise shows up in the signal may again 
impose the dynamic range limitation. We shall discuss 
these problems in the following sections. 
(b) Noise Studies of the ICP by a Fourier Transform 
spectrometer 
spectral information available from the ICP is 
encoded as a time domain signal by the FT spectrometer. 
Input radiation ie modulated at different frequencies 
according to their wavelengths. In order to obtain a 
spectrum, the time domain signal must be transformed 
through the FT algorithm. In this way, it is analgous 


to noise power spectra (83-88). Thus, one can predict 


the superimposition of the desirable signal with the noise. 


Winefordner et al. (88) have demonstrated that the 
plasma consists of 1/f, 60 Hz and a proportional noise 
centered around 200-500 Hz under the normal operating 


condition of the plasma. The exact frequency of the 


proportional noise depends on the plasma conditions. For 
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example, if the power of the plasma increases, the frequency 
of the proportional noise also increases. Work in our 
laboratory has shown similar noise behaviour of the ICP (87) 
and the proportional noise is suspected to be due to 

the asymmetrical rotation of the plasma. Because of 

the periodic nature of these noise sources, careful 
electronic filtering must be applied to the signal obtained 
from the spectrometer, or a dynamic range limitation due 

to the presence of the noise spikes will occur. The 

low frequency noises, (1/f and 60 Hz) are not of major 
concern when applying the FTI spectrometer to the ICP as they 
can be effectively filtered out electronically. The 
proportional noise spike does create a problem if one 

is not aware of its high amplitude.. An example is shown 

in Figure 42. This is a spectrum obtained using the 

Silicon diode detector. The proportional noise spike 

as indicated by the arrow can be seen near 17400 om? 
corresponding to about 400 Hz. Its presence in the 
signal will certainly impose the dynamic range problem 
for any analytical eoaiveation « The other spectrar 

lines shown in Figure 42 are due to the plasma background 
(mostly argon and oxygen lines from 400-800 nm) and 

will be discussed later in this section. Because of the 


fact that the frequency of this noise spike is fairly 


(VG) FOF 


heed dule bs meer *3 


rwoater gi: ok aang eit. 6p Sion, a att te anid 


biudagileied. Sueeiy ae st oh oa ws ordigtt al wads 


att ta bs ried 


7. on 


ant ba wd val ay sap! ora aa 
oF sgh ad ot satpdyeus- 8! eaton Senett 
tM sebsood “shines Banal te oatteton a ac ~~ f 
| nee ‘bo. won 
fod Dobe: wernoait 9 


Pe to ate 


(ise. 9a. Jabog 
i j i aaa ou 


aie! 47, 67 See 
ier gi pects v4 Be eU Vise, c w Santstt pie! th 


euv: penoty fire setlogseten | ant se 900 


tohem to. tom see) (Be 98 ‘fas aXe) ipaeton aa 


ar Bat oF notanoioses 2, adv anixsaas soci 
an 

wg? ° yt Roe hiatal the see ate 

eee irl ote £: Heys eri alge enion ano} 


oi? obey Banter 1 = SiR ein? ais 
Ph 7 a 


7 ;  . “eapfontan eboib 7 olf 
"ad OdwRE  iane mode oan verre edt yd petag ton 
sit rk Sonene mR aot. seit oou tutods ot gather geet 16) 

1 a pel 


ov sexe: Sacet’ (nda, ie ea tnéxa0 ILtiw 


aibge esto feqoad 


{exsaoge $90 eo HR 


bre hu nbaue 


107 


dOI 94} Wort paeyeutTstio estou TeuotTgzodoad jo souesezd syy SsuTZeaysn{T{TT wnayoeds *zt qyny 
f ! Id 


(,0TX) tequmusaem 
Bez2 
8t¢ S6T Sat est 


TTY rw a baa tT] 


Nay A 


Zz St 12 
(01x) epnytTTduy 


96 


Oct 


(eS Teh’ 


roqd Liew 


4 


. 


Taha 


“tT 
+ 


ws 


Sit 


77% 


j 


She aps. 


Hiskz Le 


hi 


108 


close to the frequency of the signal (1KHz-30KHz), a 
better than average electronic filter must be employed 
to prevent the noise interference on the signal. A 
variable electronic filter (Krohn-Hite, Model 3343) 
capable of 48 dB/octave rolloff per channel is employed. 
This will completely eliminate the proportional noise 
from interfering with the signal. Other means of eliminating 
the proportional noise are available. These include 
the use of a variable notch filter or. a long, extended 
ICP torch (87) in place of the conventional design. As 
mentioned before, the proportional noise does change 
in frequency as one changes the operating condition of 
the plasma. Thus, the notch filter does require frequent 
resetting if it. is being employed in the system. A 
long ICP torch simply means. an extension of the outer 
quartz tube to about 5 cm above the load coil and 
emission of the plasma passes through the quartz before 
it reaches the spectrometer. Loss of spectral signal, 
especially in the UV region is anticipated. Based on 
these facts, shnesakrohn-Hiteéseclectroniceiadtershas, been 
incorporated into the measurement electronics for the 
signal coming out from the FI spectrometer. 

As a summary, while utilizing the FT spectrometer 


as a detection system for the ICP, one has to be aware 
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of the noise present in the system as it will show up 
Superimposed with the signal spectrum. The presence of 
noise spikes in the spectrum may impose the dynamic 
range limitation on the determination of trace elements. 
The solution to this problem is the use of a good electronic 
meltemroeletiius turnytowthe problem of “dynamic range 
due to the intense plasma background and the use of a 
mixed gas) ICP asia partial solution to, this limitation. 

(c) Analytical Improvement of the ICP-FT Spectrometer 

system Using a Mixed Gas ICP 
The use of gases other than argon to sustain the 

plasma discharge can be traced back to the early stage 
Orn mhesd CPe (89-93) 409 Gases) such as helium, nitrogen, 
air and oxygen have been investigated in order to improve 
the analytical performance and to reduce the operation 
cost of the system. Recent work (92-97) has emphasized 
the use of low power mixed gas ICP's with particular 
interest in the nitrogen / argon cooled plasma. Work 
in our laboratory has shown that the nitrogen / argon 
cooled plasma is comparable to or better than the con- 
ventional pure argon plasma in terms of detection limit, 
stability, sensitivity and freedom from interference 
effects (97). Preliminary work using the mixed gas ICP 


with the FT spectrometer has demonstrated a tremendous 
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decrease in background which eases the dynamic range 
problem due to intense plasma background. In addition 
to this, the nitrogen / argon cooled plasma provides 
a more stable source when compared to the conventional ICP. 
In the following sections, the use of a mixed gas ICP 
to improve the performance of the ICP-FT spectrometer 
system will be briefly introduced. 
(i) Near Infrared Region 

As mentioned before, the detector used in this 
Pegion is the silicon diode detectoriwith a«spectram 
response ranging from 500 to 1000 nm. Within ths region, 
enisccion lines” from the alkali metals” are dominamt. With 
the use of the conventional argon ICP, plasma background 
emission (Ar, O, H lires) is very intense. Thiesrenders 
the region useless for any analytical work. The problem 
is best illustrated in Figure 43. This spectrum was 
obtained with 500 ppm sodium aspirating into the 
conventional argon plasma. As it can be seen, the 
amplitude of the sodium emission lines (589.0 nm and 
boo,.o-nm) 1s not large. In other words, the dynamic 
range problem is realized. Amplification of the signal 
is limited by the plasma background and the proportional 
noise indicated by the straight arrow in Figure 43. This spect- 


rum was obtained without the Krohn-Hite electronic i1lver and 
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this is the reason why the proportional noise spike shows 
up in the spectrum. In this situation, determination 

of sub ppm-level of the alkali metals is thus impossible. 
Introduction of nitrogen gas into the coolant stream 

of the plasma reduces the plasma background significantly. 
This is illustrated in Figure 44. This spectrum was 
obtained under the same conditions with 500 ppm Na 
aspirating. The only difference is the use of 50% 
nitrogen in the coolant stream of the plasma. In this 
case, the amplification of the sodium signal is no longer 
governed by the background. Analytical determination 

can then be performed in the sub ppm level. The multi- 
€élement alkali metal spectrum used later in this chapter 
was obtained with the aid of the nitrogen / argon plasma. 
One interesting factor that can be observed in Figure 44 
is the disappearance of the proportional noise as we 
cross over from the argon plasma to the argon / nitrogen 
Plasmas wehis further supports the fact that the argon 
plasma is rotating. The fact that the noise spike 
disappeared is likely due to the reduction in size of 

the argon / nitrogen plasma. The reduction in size 
prevents the assymmetrical rotation of the plasma and 
thus appears more stable as shown by the spectrum. In order 


to illustrate that quantitative analysis can be performed 
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with the argon / nitrogen plasma, a working curve of 
lithium covering oven a concentration range of .0.1 ppm- 
10 ppm is shown in Figure 45. 
(ii) Visible Region 

The detector uséd in this spectral region i6 a 
photomultiplier tube with a spectral response extending 
from 200-650 nm. If this detector is used in conjunction 
with the conventional argon plasma, a typical spectrum 
obtained is shown in Figure 46. Intense argon emission 
lines combined with broadband emission due to nitrogen 
species are dominant in this region. Any good analytical 
work cannot be performed due to the dynamic range limitation. 
Introduction of 50% nitrogen / 50% argon as a coolant 
gas for the ICP has a remarkable effect in reducing the 
background to a minimum. This is shown in Figure 47 
obtained using the same electronic amplification as for 
the spectrum shown in Figure 46. Note that there is 
some residual background and it is possible to further 
reduce this with the use of a 100% nitrogen cooled plasma. 
A working curve for calcium obtained using a 50% argon / 
50% nitrogen ICP is shown in Figure 48. The concentration 
of calcium ranges from 0.1-50 ppm. This concentration 
range is virtually impossible to cover with the use of 


the conventional argon ICP when coupled to an FI spectro- 
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FIGURE 45. Analytical working curve for lithium obtained 


using an argon/nitrogen ICP. 
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FIGURE 48. Analytical working curve of calcium obtained 


using a 50% argon/50% nitrogen ICP. 
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The use of nitrogen / argon mixture as the coolant 
gas for the plasma is not the only choice. An example 
of a spectrum obtained using a 50% helium / 50% argon 
PCr is shown an Figure 49. It consists of a 10 ppm 
calcium signal and a number of intense hydrogen emission 
lines. The presence of the hydrogen lines can be due 
to the hydrogen gas impurity in the helium gas supply 
or can be simply due to the hydrogen content of the water. 
In either case, it can be easily removed by filtration 
or desolvation techniques. In this way, a super clean 
plasma background can be obtained. The only major 
disadvantage of the use of a helium / argon ICP is the 
Mastaer lity of the plasma. it%is illustrated here 
simply to demonstrate that the argon / nitrogen ICP may not be 
the ultimate choice for the FI spectrometer system; 
further investigation is necessary. Besides, the optimum 
observation height, gas flow rates, power levels have 
still to be determined for the nitrogen / argon ICP. 
This investigation is best carried out using a photodiode 
array spectrometer and is currently underway in our 
laboratory. The work presented ‘pera by no means represents 
the optimum signals that can be obtained by the system. 


Further investigation must be carried out to fully extend 
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the capability of the ICP-FT spectrometer system. 
(iii) Ultraviolet Region 

Generally speaking, as pointed out by Hall (33), the 
conventional argon ICP provides an exceptionally clean 
Spectral region for analysis with a solar blind PMT as 
a detector. The use of argon / nitrogen ICP will 
increase the molecular band emission due to nitrogen 
Species. Thus, normally a: conventional argon ICP is 
employed. A working curve of magnesium with concentration 
ranging from 0,1-50 ppm is shown in Figure 50 for 
a) lustration. 

(ad) SConclusion 

The feasibility of the ICP-FT spectrometer system 
for good analytical work is hindered by the dynamic 
range problem present due to the multiplex property of 
the spectrometer. Through the use of a mixed gas ICP 
such as a 50% nitrogen / 50% argon plasma, the dynamic 
range limitation due to the strong background emission 
of the plasma has been reduced to a minimal level. 
Depending on the spectral region we are interested in, 
different plasma excitation sources can be utilized to 
reduce the background. This is summarized in Table 4. 
In addition, the argon / nitrogen ICP is a more stable 


Source when compared to the conventional argon ICP. 
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FIGURE 50. Analytical working curve of magnesium obtained 


using a conventional argon ICP. 
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However, other means of reducing the plasma background 
Such as the use of a pure nitrogen plasma proposed 

by Barnes and Meyer (98) or the use of a microwave 

induced plasma which has a fairly clean spectral background 
(99) as compared to the conventional argon ICP cannot 

be overlooked. But, this will not improve the dynamic 
range problem imposed by the matrix of the sample. The 
peeeoi,Opeical filters combined with the availability of 

a good selection of spectral lines due to the wide 
spectral coverage of the spectrometer can in some way 
overcome the problem. For very complex matrices, 

Optical isolation of the strong emission lines from 

weak emission lines must be performed. Techniques such 

as the selectively modulated interferometric dispersive 
spectrometer (100-102) or dual beam interferometer (103- 
105). can, pe. employed... In the, former, case, dhe tixed 
mirror of the interferometer is replaced by a diffraction 
grating so that only a narrow band of spectral information 
is modulated and utilized by the. spectrometer. In this 
case, the dynamic range problem will be completely 
eliminated. However, the multielement capability of the 
system is reduced or completely lost. However, the precise 
wavelength calibration inherited in the system can still be 


fully utilized. In other words, a precise slew scan type 
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Spectrometer can be obtained. 

In the case of a dual beam interferometer, the 
complementary outputs of the interferometer are used 
to cancel out the very intense emission signal present 
in the system. Thus, the weak emission spectral infor- 
Mation can be detected. However, practical difficulties 
in exactly cancelling out the strong emission signals might 
arise. 

In summary, the application of the ICP-FT spectro- 
meter system to real sample analysis cannot be overlooked 
and the system will certainly extend to its full potential 


iaethesdynamic ranger dimitation can! berfullyoovercome. 


D. Cross-correlation Masks for Processing Interferograms 
The same concept of correlation techniques as 
applied to the spectral signals obtained using a photo- 

diode array spectrometer can be utilized in a similar 
fashion for data obtained using a Fourier transform 
spectrometer. However, if correlation analysis is to 

be performed on signals in the spectral domain, at least 
two FFT's will be involved to convert the inter- 
ferograms into spectra, one for the experimentally 
generated cross-correlation mask and one for the 


experimental signal. This is certainly not a very 
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efficient way of performing correlation analysis. There 
is no particular reason why cross-correlation techniques 
cannot be applied to time domain signals even though humans 
tind these signals difficult to directly interpret. Hall 
(33) reported the use of experimentally generated 
interferograms as cross-correlation masks for extraction 
of analytical information from signal interferograms. In 
tiie case, ‘Selection of specific spectral features in 

the cross-correlation masks cannot be accomplished 

eacily, in Contrast’ to the correlation masks used in 
spectral domain, by means of threshold levels. Threshold 
levels can be implemented in a similar way only if 

the cross-correlation masks are Fourier transformed, 
threshold imposed and inverse Fourier transformed back 

Se thmegdomainesignals; winhoaddi tipni nidrtsrveryedifficult 


to obtain a very "clean" cross-correlation mask containing 


only the specific spectral features required experimentally. 


In most cases, background emission features are also 
included in the cross-correlation masks. Elimination of 
overlapping lines is also difficult. Further, phase 
differences between the cross-correlation masks and 
signals were corrected by performing a partial full 
cross-correlation function. Relative displacement 


(typically, ten data points along positive and negative t 
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displacements) between the signals and the correlation 
masks was performed in order to obtain the desirable 
correlation data. This increased the computation time 
for the cross-correlation operation. In this chapter, 
software generated cross-correlation masks will be 
utilized. This is possible only because of the precise 
wavenumber axis provided by the Fourier transform spectro- 
meter. 

Let us recall that the cross-correlation function 
at the t = 0 point can be expressed mathematically as 


the following: 
Gigs > a0, ie) 


where C is the cross-correlation data point at t = 0, b 
is the digitized signal which is an interferogram in 
ters case and a is the” cross-correlation mask. The 
correlation data point is obtained by the summation of the 
products of the point-to-point multiplication between 
the signal and the mask and it will contain the desirable 
analytical information as selected by the mask. 

If monochromatic radiation is directed into a 
Michelson interferometer, a cosine waveform interferogram 
results. An example using a He-Ne laser as the source 


of radiation is depicted in Figure 51 with the dotted 
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He-Ne Laser 
(632.8 nm) 


Michelson 


interferometer 


FIGURE 51. Interferogram resulting from a He-Ne laser. 
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line representing the zero path difference between the 

two mirrors in the interferometer. In an ideal case, that 
is no phase shift is present, the cosine wave originates 
at the zero path difference point. When the plasma 
system is coupled to the interferometer for simultaneous 
multielement analysis, the interferogram obtained is 
simply a waveform resulting from the summation of 
different cosine waves, each cosine wave corresponding 

to a particular spectral line. With a laser referenced 
interferometer, the exact waveform of the Signals can 

hemp caicvcd, that is to say,.for every spectral line, a 
precise cosine wave can be generated in the computer 

to be used as fa cross-correlation mask. The idea of/ cross- 
correlation techniques as applied to interferograms is 
illustrated in Figure 52. The cross-correlation operation 
between an interferogram and a cross-correlation mask 

is simplified with the use of cosine waves (i.e. mono- 
ehromatic svenais) . When correlation is performed 
between a signal and a mask that exactly match with each 
other, as shown in Figure 52(a), the product of every 

pair of data points will always turn out to be positive. 

A high correlation data point will result indicating 

the similarities between the signal and mask. As shown 


in Figure 52(b), if the same correlation mask is used for 
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(a) 


(b) 


FIGURE 52. Cross-correlation between cosine waves of 


(a) same frequency and (b) different frequencies. 
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signals that do not correspond to the same frequency, 


the summation of the point-to-point products will average 


CUtebOrd Very; LOW correlation data point. ‘Thus, cross- 
correlation analysis is analogous to the classical 
Fourier transform method, with only limited desirable 
spectral information being processed. Hall (33) has 
made an excellent comparison between correlation and 
Fourier transform processing methods and his work 
should be consulted for more details. 

As can be seen from the figure, the phase relation- 
Ship between the cross-correlation mask and the signal 
is very important in order to get the right information 
euu jot the sienal:, A phase correction routine is 
incorporated into the data processing step. It simply 
involves the generation of the appropriate cosine and 
Sine waves as two different cross-correlation masks. 
Correlation analysis.1e then penformed between these 
masks and the signal individually. Two correlation 


data points will be obtained, one for sine and one for 


cosine. The root mean square value of these two correlation 


data points will be the phase corrected correlation data 
point +») The; cosine: and ‘sine cross-correlation masks 
which are essentially synthetic interferograms are 


generated by the following expressions: 
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where v is the wavenumber of the spectral line in interest, 
4S is the sampling interval and n is the point number 
of the interferograms. In our interferometer system, 
a He-Ne laser was employed to provide a sampling interval 
or 0.63281. 

A double sided interferogram is normally required 
for phase correction purposes when the interferogram 
me processed through the FFP route. This is not necessary 
in the case of cross-correlation techniques. Phase 
correction is accomplished by the use of cosine and sine 
waveforms. Thus, a single sided interferogram is 
Suriicient. This immediately reduces the number of data 
points to be processed by correlation analysis in half 
as compared to the FFT route while maintaining the same 
degree of resolution. In this study, interferograms after 
the zero path difference point during the forward mirror 
movement will be utilized as opposed to the double sided 
interferograms. 

In addition to the cosine and sine cross-correlation 
masks, the use of square wave synthetic interferograms 


as cross-correlation masks is also investigated. Square 
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synthetic interferograms are very much the same as the 
cosine and sine synthetic interferograms except they 

contain values of 1,0 and -1 only. They are obtained 

by converting all values of the cosine or sine synthetic 
interferograms that are greater than zero to 1's 

and -1 for values that are less than zero. That is 

to say, the amplitude information is discarded while 
retaining the frequency information of the cross-correlation 
masks. By doing this, the correlation process is simplified, 
as Multiplicatroneissreplaced by.addition or subtraction. 
This will certainly speed up the correlation process 

and in addition, less memory space is required to store 
these cross-correlation masks in a computer. In the 
following sections, the cross-correlation technique 

using both the analog (cosine and sine) and square 

synthetic interferograms as cross-correlation masks will 


be described. 


E. Elemental Analysis Using Cross-correlation Techniques 
Qualitative Aspects: A multielement interferogram 

containing sodium, potassium, rubidium, cesium and 

lithium information is shown in Figure 53. Obviously, 


the interferogram has no meaning at all as seen by 


human eyes. This interferogram was Fourier transformed 
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Haro its spectral domain as shown in Figure 54. In the 

spectral domain, it would then be possible to tell 

which elements are present. In this section, cross-. 

correlation techniques for automatic elemental analysis 

of the signal in the Fourier domain will be illustrated. 
As the multielement alkali metal interferogram 

was stored in the computer memory, synthetic interfero- 

grams or correlation masks, each corresponding to a particular 

wavelength of an element, were generated and correlated 

with the signal. The correlation masks of seventy 

elements were generated using the most sénsitive ICP 

lines as reported by Boumans (106). These lines are 

included in Table 5. Each element has been assigned a 

code number arranged alphabetically for simplicity. There 

is no reason why the atomic number of the elements 

cannot be used. The present coding system was adopted 

simply because the ICP spectral lines of each element 

are listed alphabetically in Boumans paper. It must 

be kept in mind that the most sensitive ICP lines presented 

in Table 5 may not be true for the ICP-Fourier transform 

spectrometer system, especially when a mixed gas ICP 

is utilized. However, the table should be good enough 

for preliminary work. For future consideration, an 


atlas of ICP lines based on the interferometer system 
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TABLE 5. Most sensitive ICP emission line for 70 elements. 


Element Code Number Wavelength (nm) Element Code Number Wavelength (nm) 
Ag 1 328.068 Na 36 588.995 
Al 2 396.152 Nb 37 202.932 
As 3 197.197 Nd 38 430.358 
Au 4 201.200 Ni 39 231.604 
B 5 208.959 Os 40 206.721 
Ba 6 455.403 P 41 PLS) lente! 
Be 7 234.861 Pb 42 405.783 
Bi 8 195.389 Pd 43 340.45€ 
e 9 247.857 Pr uy b17.93€ 
ap 10 393.366 Pt Ls 204 .937 
of} ai 228.802 Rb 46 760.023 
Ce a2 395.254 Re 4? 227.525 
Co 13 228.616 Rh 48 343.489 
Cr 14 205.552 Ru 4g 240.272 
Cs 15 852.124 Sb 50 206 .$33 
Cu 16 324.754 Se 51 361.364 
Dy 17 353.170 Se 52 196.026 
Er 18 390.631 csi 53 251.611 
Eu 19 420.505 Sm 54 356.627 
Fe 20 259.940 Sn 55 197.080 
Ga 21 294.364 Sz 56 407.771 
Gd 22 342.247 Ta 57 214.687 
Ge 23 199.824 Tb 58 350.917 
Hf 24 202.818 Te 59 200.202 
Hg 25 253.652 Th 60 401.913 
Ho 26 345.600 Ty 61 334.941 
In 27, 325.609 TR 62 351.924 
Tr 28 216.942 Tm 63 346.220 
K 29 766.490 U 64 385.958 
La 30 394.910 V ; 65 309.311 
Li 31 670.784 W 66 202.998 
Lu 32 261.542 Y 67 371.030 
Mg oD 279.553 Yb 68 369.419 
Mn 34 257.610 zn 69 206.200 


Mo 35 202.030 he 20 339.198 
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Similar to the one obtained by Fassel et al. for dis- 
persive spectrometers (107) should be obtained. 

The results of the cross-correlation process for 
all seventy elements are plotted in Figure Bee This 
is applot of the cross-correlation functi6én at t = 0 
point versus an element axis where the elements are 
represented by their respective code numbers. From the plot, 
it can be seen that all the elements that were present in 
the sample show relatively high correlation results as 
compared with the others. In addition, indium, which was 
absent from the sample, also shows relatively high correlation 
datas This is because the indium line (325.61 %nm), if 
present, will overlap with the lithium line (670.78 nm) 
because of_alaising (32). However, the presence of indium in 
the sample is unlikely because the response of the silicon 
diode detector used to obtain the signal is minimal near 
300 nm. Similar results were obtained using square synthetic 
interferograms as the cross-correlation masks. This is 
shown in Figure 56. The selectivity of the analog and the 
square synthetic interferograms appears to be comparable. 
The only difference is the slight increase in the absolute 
magnitude of the correlation data obtained using the square 
synthetic interferograms as the cross-correlation masks. 


This is because the square synthetic interferograms do 
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not contain any value between O and -1; O and +l. The 
absolute magnitude of the correlation results represents 
the quantitative information about the sample analysed. 
whis will be further discussed later in this chapter. In 
this way, the relative magnitude of the correlation 

data can be utilized to indicate the qualitative 
information of the signal. It can be implemented in 

an automatic fashion by a computer for spectral analysis. 
This is analogous to spectral search methods used in 

the infrared spectroscopy for compound identification. 

A spectral match between the signal and a library of 
spectra is performed in the spectral domain. With the 
correlation technique, spectral analysis can be accomplished 
in the Fourier domain utilizing a library containing 
various cross-correlation masks for each element. 

The correlation technique was further investigated 
with the use of a more complex spectrum. The signal 
containing information about some transition metals 
is shown in Figure 57. These transition metals were 
cobalt, iron, nickel and vanadium. For the sake of 
illustration, this interferogram was Fourier transformed 
into the spectral domain and is shown in Figure 58. A 
fairly complex spectrum was obtained. The pure spectra 


of these elements, are shown in Figure 59 for reference. 
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Note that in the case of correlation analysis neither the exact 
position of the spectral lines nor the degree of spectral 
overlap need be known. Cross-correlation analysis was 
applied to the time domain signal using cross-correlation 
masks generated according to the wavelengths listed in 
Table 5. These results were plotted as shown in Figure 
60. The information obtained is erratic because neither 
sodium nor palladium was present. Besides, two out 

of the four elements that were present in the sample 

do not show relatively high correlation results as 
compared to the others. This problem was caused by the 
poor selection of wavelengths used to generate the cross- 
correlation masks. With the use of«the=solar blind 

filter and the visible PMT combination as the”detector 

am obtaining this signal, the spectrateresponse of the 
system was limited to within 240 nm to 360 nm region. 

The correlation analysis was then repeated using the 
cross-correlation masks generated by the wavelengths 
listed in Table 6. This table contains the most sensitive 
line for each element within the spectral window of 240 nm 
to 360 nm. Again, the wavelengths were extracted from 

the table published by Boumans (106). For some elements, 
there are no reasonably sensitive spectral lines within 


this region, the spectral line that is closest to the 
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TABLE 6. Most sensitive ICP emission line for 70 elements 


Element 


Ag 
Al 


within the spectral window of 240-360 nm, 
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Wavelength (nm) 


328.068 
309.271 
278 020 
267.595 
249.773 
234.758 
313.042 
306.772 
247.857 
317.933 
3L6 620 
394.275 
345.350 
283.563 
455-531 
324.754 
3534170 
337.271 
381.967 
259.940 
294.364 
342.247 
265.118 
339.980 
253.652 
345.600 
325.609 
254.397 
404.414 
394.910 
323.263 
261.542 
279 +553 
257.610 
313.259 


Element 


Na 
Nb 
Nd 


Code Number 


36 
37 
38 
39 
4O 
41 
42 
43 
Qu 
45 
Lé 
47 
48 
4g 
50 
51 
52 
53 
54 
55 
56 
5? 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 


Wavelength (nm) 


330. 
-418 
+333 
-166 


309 
386 
310 


253. 
253. 
»199 
458 
390. 
G45 
420. 


280 
340 
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region was chosen to complete the list. The results are 
Shown in Figure 61. All four elements that were present 
in the sample show relatively high correlation, in 
addition, aluminium and palladium also show strong 
indication of their presence. This is because the 
aluminium line at 309.27 nm overlaps with the vanadium 
trnevay 309.31 nm and in the case of pailitadium (340.46 nm), 
there is a spectral overlap with a cobalt emission line 
Berot0Oa5bnm. It must a@lso@be pointed™out=thatethe 
relatively high "background" for the rest of the elements 
tested is mainly caused by the complexity of the signal. 
Due to the large amount of emission lines that are 
available from the four elements present in the sample, 
spectral overlap is severe. 

In order to. distinguish the elements. that were 
present in the sample from the elements that were not, 
but showed high correlation, a secondary correlation 
process is necessary. This simply means the utilization 
of more spectral lines of the suspected elements for 
correlation, analysis @elnethis=case two additional 
spectral linesMof each suspected element were used for 
' the generation of the additional cross-correlation values. 
There is no particular reason why a total of three 


spectral lines were chosen. In complex situations where 
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spectral overlap is more severe, more than three spectral 
lines for each element can be used to generate the cross- 
correlation masks for analysis. The wavelengths of 
three spectral lines for each element are shown in Table 
7» Each spectral line has been assigned a code number 
as indicated in Table 7. Each of these spectral 

lines was used to generate a cross-correlation mask for 
analysis and thus each element has three correlation 
data points corresponding to each spectral line. The 
results of these are plotted as shown in Figure 62. In 
the case of aluminium and palladium, only the spectral 
lines (aluminium 309.27 nm and palladium 340.46 nm) 

that overlap with the lines originating from other elements 
which were present have high correlation data points. 

For the other spectral information tested for these 
elements none of them yields a high correlation result. 
In contrast, all four elements that were present in the 
sample show relatively high correlation for all the 
spectral information tested. In this way, one is-able 
to pinpoint exactly which elements are present in the 
Signal and automatically identify spectral overlap 
problems. As mentioned above, the absolute magnitude of 
the correlation data indicates the amplitude of the 


particular spectral information in the signal. 


igegn sae arr amet EO ease so at « 
=gRoun. sit avirtongyy: ‘ath wea od an soamante 
to atye lever Pe weit <Or'e 
ehdet ont twas ais: — fos 2% mere 
tou a gg Hi calli ceaid Hat ‘sokt Laat 
erase Sesths a: doed-/ ( ofdel aks 
dt wei: nb Pe ekesio eee m, ae1siey of bad 
| Hichyon sere ashi! ‘3eomete goes euAs bes 
eat. a . - Lasrtoagae tiakap OR ‘grhamigeer toe ate 
¥E; (SS et yh i Rohe, be nadesta ate sasdyt bad ark 
cettoven add vido JAD ERa Dey one avietonla to" 2880 
nyt BE. debe tieg es tie FS. WOE a Ered tute} 
ornatela cotta Moet galianitgien gomdt: ‘oat sorb galt: 
ateiod tee neiteieriwe. Lai wea tneasiq. § 
saat “alt Doitee, oleihpaaieanes Niewtosqa vedio 
tes OstHtbetabO: ‘Aiilee = Sbbb Te Heds ‘to. sor 
sdk 2 sane stongaay nteumete Hane Oem a in 
| ay Cie ‘ph soho. statis aaah 
pide 8 ane ond Sannin 
BAS ast tress | : anaae 
_ gatisve mati, ea vee obi pte ponge b x 
1? edu Enger orukaeete! iit avods +o iat 8A ‘ameLd 7; 
sit ‘to obit Lai oie eal tab tteaiersey on 


Kenge itt — becoaaae “aluote 


, a 


isms 
Seok 


P51 


TABLE 7 -Some-1CP+emission~Lines. for—aluminium,.cobalt, iron, 


nickel, palladium and vanadium (240 nm - 360 nm). 


spectral Line 
Element 
Code Number 


Wavelength (nm) 
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FIGURE 62. Secondary cross-correlation analysis results 
for (a) aluminium, (6) cobalt, (c)eiron, 
(d) nickel, (e) palladium and (f) vanadium 


utilizing analog cross-correlation masks. 
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An alternative of performing the secondary correlation 
analysis is to combine the correlation results obtained 
from more than one correlation mask for each element. The 
results are shown in Figure 63. In this case, a code 
mumber Of}; one represents the correlation data point 
obtained from the correlation mask generated using spectral 
line #1 for that particular element. A code number of 
two represents the summation of the correlation data 
points obtained from the correlation mask generated using 
spectral lines #1 and #2 for that particular element and 
SO On. 

As shown in Figure 63, in the case that high 
correlation is due to spectral overlap, the summation of 
all the correlation results obtained from three independent 
correlation analysis remains fairly constant as in the 
case of aluminium and palladium. Cobalt, iron, vanadium 
and nickel all show a staircase fashion increase in the 
correlation results. Thus, this can also be utilized as a 
mean for qualitative analysis. All the above analysis 
have been performed using the analog cross-correlation 
masks (i.e. cosines and sines). Similar results using 
the square synthetic interferograms as cross-correlation 


masks can be obtained. The results are illustrated in 


Figure 64,65 and 66. 
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FIGURE 63. 
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FIGURE 65. Secondary cross-correlation analysis results for 
(Calealuminiumy (b) Cobaltem(e)udarony) (d)_nickel, 
(e) palladium and (f) vanadium utilizing square 


cross-correlation masks. 
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In situations where spectral analysis is difficult, 
cross-correlation analysis can be very helpful. Most 
of all, correlation analysis can be implemented essentially 
in an automatic fashion in a small computer system. Qual- 
itative information can be obtained within a short 
period of time by performing the analysis using cross- 
correlation masks stored in a computer. In this case, neither 
knowledge of the spectral axis nor the degree of spectral 
overlap need be known once the library of/correlation 
masks has been established . 

Quantitative Aspects: As mentioned before, the 
absolute magnitude of a correlation data point represents 
the quantitative information at a particular spectral 
wavelength. Linear working curves for lithium with 
concentration ranging from 0.1 ppm to 10 ppm are shown 
in Figure 67. They were obtained using the analog and 
square synthetic interferograms corresponding to the 
dithium line au 670.78 nm’as correlation masks... As 
indicated by the slopes of the graph, the results 
obtained using square synthetic interferograms as 
correlation masks have a better sensitivity over those 
obtained using analog synthetic interferograms. Both 
methods do demonstrate their abilities to extract 


quantitative information from analytical signals. 
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FIGURE 67. Linear analytical curves for lithium obtained by 


(a) .square and (b) analog cross-correlation masks. 
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nie elrcet (of combining correlation data obtamned 
from different correlation masks, each corresponding 
to the quantitative information at a particular spectral 
Wavebensth, {is illustrated in Figure 68. This is a plot 
of conrelation results versus the concentration of 
fwekes With|concentration ranging from 0 t6v500 ppm. As 


moresspectral information is utilized, the sensitivity 


as indicated by the slopes of the graph increases. However, 


One Gannot Simply utilize all the spectral information 
available in the signal as spectral interferences may 

be present. This is best illustrated with an example, 
Getermmination of vanadium and nickel. Among the more 
sensitive lines of vanadium®and nickel, the vanadium 
Pine-at 310223 nm overlaps with the nickel line at 310.17 
nm. If the three most sensitive spectral lines of 
Vonaqdium andg=nickev are ulimized’ to establish an analytical 
Calibration| curve, the results illustrated in Figure 69 
Will be obtained. Correlation results using nickel and 
vanadium cross-correlation masks on solutions containing 
a constant amount of vanadium and an increasing amount 

of nickel are plotted as shown in Figure 69(a) and 69(b). 
Obviously, the vanadium cross-correlation masks are 
extracting nickel information because of the spectral 


overlap between the two elements. Thus, one has to be 
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Linear analytical curves for nickel obtained by 
cross-correlation analysis with analytical 
information corresponding to (a) 310.17, 341.48, 
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FIGURE 69. Detection of nickel and vanadium in the presence 
of each other with cross-correlation masks 
corresponding to (a) 310.17; G4T.48, 305.08; 

Cp) B09 sei ol 07, SiO, 23 2nm-. 
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careful about the choice of correlation masks. Upon elimination 
of the 310 nm spectral line for both nickel and vanadiun, 
more desirable results are obtained as shown in Figure 70. 
Correlation analysis on vanadium shows a constant level 

of correlation data. In the case of nickel, a linear 
calibration curve with zero intercept is obtained. If 

a plot of correlation data obtained for a particular 
wavelength against the concentration shows a non-zero 
intercept, spectral interferences should be suspected. 

A final example is given for the determination of vanadium 
in the presence of cobalt, nickel and iron. A linear 
calibration curve for vanadium concentration ranging 

frome vo 10-ppm is obtayied as illustrated in Figure 71. 
The correlation data correspond to the vanadium information 
an 309.3. nm. 

Precision is also an important consideration when 
dealing with the correlation analysis. A precision study 
On a ysolutiom ot 10 ppm vanadium was carried out, § ihe 
data arcUlisted=inyrapte-O>—Thesernumbersrare tie percent 
relative standard deviation of the 10 ppm vanadium signal 
over ten readings. The conventional peak analysis and 
correlation analysis based on both the analog and square 


correlation masks give comparable results. 
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FIGURE 70. Detection of nickel and vanadium in the presence 
of each other with cross-correlation masks 
corresponding to (a) 341.48, 305.08; (b) 309.31, 


SL Ye07 “nin. 


it obe eon ce roeg 
k a , ; : <\n 


Pee = n Les - 7 
Sovessed Smt rk. meLpal 


Lea: 


meme 


7) - : i i my if - 


7 Bibel = - Ms ah 
2teem noitsis woereaos 
$c. 08 (a) 780.2007 Gee 


a 


{ 


a 
oo ie 


165 
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FIGURE 71. Analytical working curve for vanadium obtained 
in the presence of cobalt, nickel and iron by 


cross-correlation analysis. 
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F. Conclusions 

A final comparison between correlation analysis and the 
Fourier transform data processing method as applied to 
the ICP-FT spectrometer system is depicted in Figure 72. 
All the relevant analytical information can be obtained 
by correlation analysis without the necessity to 
convert the spectral data into the human interpretable 
spectral domain. Moreover, correlation based data 
processing of interferometric signals is simple when 
compared to the conventional fast Fourier transform 
data processing method and yet is a powerful technique. 
Because of its simplicity, the need for a large computer 
to carry out data processing for long interferograms 
disappears. Furthermore, the correlation analysis 
can be implemented essentially in an automatic fashion 
by a computer to perform qualitative and quantitative 
analysis within a very short period of time. The 
flexibility in the generation of cross-correlation masks 
fully utilizes the wide spectral coverage of the ICP- 
Fourier transform spectrometer system. Different spectral 
information can be obtained simply by using different 
cross-correlation masks. Background correction by 
correlation analysis is also possible. This can be accom- 


plished by using an appropriate correlation mask that 
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resembles closely in frequency the correlation mask 
used to extract the desirable analytical information. 

With the advent of electronic multiplier chips, 
correlation analysis which is essentially a process of 
multiplications and accumulation can be accomplished in 
real time. In the next chapter, a correlation based real 
time data processing system using a microcomputer as a 
controller designed for the ICP-Fourier transform spectro- 


meter system will be described. 
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CHAPTER V 


A Real Time Correlation Based 


Data Processing System for Interferometric Signals 


A. Introduction to Real Time Data Processing Systems for 
Interferometric Signals 

Recent hardware developments have resulted in the 
appearance of devices capable of essentially real time 
Fourier transformation of electrical (66) and optical (67) 
Signals. The electro-optical system is particularly interest- 
ing in that a two dimensional Fourier transform of an optical 
image can be generated in a matter of milliseconds. High 
speed fast Fourier transformation is also available as 
computer firmware (68). Though the number of data points that 
the firmware can handle is limited at the present stage, the 
ultimate goal of performing fast Fourier transformation for 
a large block of data essentially in real time is not far 
away. These and future developments are likely to revolution- 
ize the present data processing methods involving Fourier 
+ranshormatlontmApplacationstofSstheserdevices tomrealrtime 
data processing for the ICP-Fourier transform spectrometer 
system may be possible in the near future. In this chapter, 
a real time data processing system based on the correlation 
approach will be described. It employs the use of the cross- 
correlation function at t = 0 point, in consequence of the 


success of the software implementation of the correlation 
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technique upon stored interfergrams described in Chapter IV. 
As the fast Fourier transform algorithm for processing 
interferograms is not involved when correlation techniques 
ere employed, the useiof a microcomputer 1S suiticient for 
the correlation mathematical computation, even in the case 
of long interferograms. This can be achieved either in a 
software or hardware approach. In a software approach, 
interferograms stored in computer memories can be processed 


in exactly the same fashion as described earlier. However, 


correlation which is essentially a process of multiplications 


followed by accumulation can be hardware implemented more 
efficiently by utilizing high speed multiplier - 
accumulator electronic chips. The first prototype of a 
microcomputer based data processing unit designed for the 
ICP-Fourier transform spectrometer system will be presented. 
The basic idea is to use a microcomputer as a controller 
with a multiplier - accumulator electronic chip acting as 
the actual data processor. Cross-correlation masks are 
stored in computer memories and sent out to the multiplier - 
accumulator at appropriate times for processing. Multiplica- 
tion between the acquired signal and the correlation mask 
followed by accumulation is performed in real time. Due to 
the fact that phase correction is necessary, two individual 


scans of the same signal are required for one correlation 
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data point (one for the cosine multiplication and one for 
the sine multiplication) because only one accumulator is 
available in the multiplier-accumulator chip. 

This system is not meant to be a definitive configurat- 
ion but wassonly-.constructed to test the basic idea of a 


Simple hardware correlation processor. 


B. Rockwell AIM 65/TRW Multiplier-Accumulator Data Processing 
system 

A block diagram of the experimental system is shown in 
Presume 73. Lt Can be subdivided into the following: 

(a) Rockwell AIM 65 microcomputer. 

(bo) Interferometer»control and measurement electronics. 

(ce) Star pulse ands elock control circuitry. 

(d) Analog to digiwer eonversion circuitry. 

(e) TRW multiplier-accumulator. 

(f) Data multiplkexine circuitry. 

The microcomputermeacts asva controller for she wwhoie 
system and as the source of the correlation masks necessary 
for the data processing. The analog signal (interferogram) 
from the interferometer system is digitized by the analog 
to digital conversion circuitry. The computer keeps track 
of the end of conversion pulse for every signal and as 


Soon as the analog to digital conversion is finished, the 


= 6 Oe ce 
a 5 i: Fi 


vot sme -bets ao Peer: ra 
ai yotelumutan: a eo 


wlke ‘ - 4 > we 
Woiis + c nts y= 
; Ly 


aitivol int aAr ost cance ao a  Gaaen* va | 


‘oats Jnewstifeaser saa’ igs J 


y i - wy 
a Lot ‘ afi? hat ei noabusiiann vw nee one 5 \ au 
f is - 


vineadoe Size Saireahols — ‘ i 
(me2ao'es Retin: } iain sf a ait 
abides add vid siia bbe é 
doas? syead 12 tuqnod sho) ga 

aa Bae facyts eee z oes 


Su? .oetainét a2 axpizeeer 


173 


‘uayshs Sutsseooad eyep 


lozeTnunsoe-JatTdtytTnw myL/teyndwosozotw $9 WIV TTeMyoOY ey JO weisetp yootqd °*€24 AYNDId 


Bye - 
Auynoutg 1€Q 14-8 


Bujxedinjinw nil 


49019 
| _$}044u09 | j2y161g sese7 
sQB|NWNDDy EE syNdwos0s9\~Wy Aayinoug 3 Baa 
/aNdnINW SAS|Nd YO0|D SO WIV [013U0D 49019 JEUDIS JU04}98/3 


|jEMxDOH SNOUDIYOUAS| JOJOWOIBP9}Uj 


PUBWWOD | pue easing BIS 


HSE 
U0}}8|90D 


yeuBis 14617 SUM 


UO|SiB@AUOD 
jo pug 


Auynds1D 
UuOISJBAUOD 
j2316)q 0) Bojeuy 


puewwoy 


Byeq wa-Ol Tenuo5 


jeubis Bbojeuy 


a _— - = ~~ é tein 7 Ce. 
S pete ee ee «6 eee 
a : ay wae a ae) ae ee ——-- 


i) WiA isewsoot ety “to matgekhcleokt «Ey 2arprt 


, x 


wataes arivasnoty = 


appropriate correlation mask data point is fed into the 
multiplier followed by the appropriate control and clock 
Signals. At the end of a scan, as recognised by the 
computer, control signals are sent to the multiplexing 
circuitry in order to multiplex the final data into the 
computer in 8-bit data segments. This result will then be 
manipulated as required and output to display. 

vert spulsejand=acieckicontraovocincurtry which,is 
controlled by the signals coming from the interferometer 


electronics is incorporated into the system so that the 


data acquisition/processing will start at the right moment. 


Each of these subsystems will now be described in 
more detail in the following sections. 

Rockwell AIM 65 Microcomputer: The "heart" of the 
data processing system is a Rockwell R6500 Advanced 
Interactive Microcomputer (AIM 65). It consists of two 
modules; a master module with a 20-column thermal printer, 
a 20-character display, and the microcomputer components; 
a standard keyboard module for communication between the 


computer and the operator. The R6502 Central Processing 


Unit (CPU) of the microcomputer operates at 1 MHz. Software 


support of the AIM 65 computer consists of an 8K MONITOR ~ 
which includes an EDITOR; 8K BASIC and 4K ASSEMBLER. A 


typical computer like this with 4K RAM (Random Access 
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Memory) costs less than $1000 and thus Significantly reduces 
the initial cost of the data processing system. 

The interface of the computer with the rest of the 
system was accomplished through the R6522 Versatile 
Interface Adapter (VIA) which is available to the user at 
the Application Connector of the computer. The VIA consists 
of the following: two 8-bit I/0 ports (A and B), four 
peripheral control/status lines (two for each port), two 
16-bit counters/timers and an 8-bit shift register. Only the 
I/O ports and the control lines were utilized. The detailed 
configuration of the VIA is shown in Figure 74. Port A I/0 
lines (PAO-PA7) were used for data input into the computer. 
The accumulator output from the multiplier-accumulator chip 
was multiplexed into the computer in 8-bit segments through 
Port A. There are two peripheral control/status lines at 
Port A, namely CAi and CA2. The CA1 line was set up so 
that the end of conversion signal (1-0 transition) from the 
analog-to-digital converter could be monitored. A positive 
pulse which acted as a multiply command signal to the 
multiplier-accumulator was generated at the CAz2 line. Port 
B I/O lines (PBO-PB7) were used for output control signals 
and correlation masks. Correlation masks were sent from 
the computer to the multiplier-accumulator through the PBO 


and PB1t lines. Accumulation and output register tristate 
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controls of the multiplier-accumulator were set by the PB2 
and PB5 lines respectively. For multiplexing accumulator 
data into the computer, two control lines for the 
multiplexing circuitry were necessary. These were controlled 
by the computer utilizing the PB7 and PB6 1/0 lines. In 
order to have the whole system running, a CLEAR command 
followed by a GO command which were software controllable 
must be issued through the PB3 and PB4Y lines respectively. 
Only one of the two control/status lines at Port B were 
utilized. The CB2 line was set up to generate a positive 
pulse as a data loading command signal for the multiplier- 
accumulator. 

The exact sequence of events during the operation of 
the system is shown schematically in Figure 75 and should 
be consulted durimge the discussion,of the rest) of the 
system. 

Interferometer Control and Measurement Electronics: The 
interferometer electronics can be subdivided into two parts: 
a controlling part and a measuring part. The control elect- 
ronics have been described in detail by Yuen (32) and later 
py Hall (33). Three digital signals from the control 
electronics were utilized. They were the white 
light pulse, the synchronous signal and the laser digital 


Clock. Because of the fact that a single sided interferogram 
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is sufficient for correlation analysis, the control 
electronics can be substantially simplified. The phase lock 
loop used to stabilize the mirror drive, the white light, 
the synchronous signal and the laser signal are sufficient 
for the full operation of the correlation data processing 
system. However, the three signal outputs were obtained 
directly from the control electronics circuitry to 
maintain the proper operation of the original system. 

The white light pulse, together with the synchronous 
Signal and the laser digital clock signal were connected 
airectlyninto)thecstart pulse.and.clock, control circuitry. 
The white light pulse was used as the start pulse for the 
data acquisition and processing. The interferogram obtained 
after the zero path difference point between the two mirrors 
in the interfermeter as indicated by the appearance of 
white light pulses during the forward movement of the mirror 
was used as the signal for processing. This was recognized 
by the electronics with the aid of the synchronous signal. 
The laser digital signal was used as the conversion clock 
for the analog to digital conversion circuitry. 

The measurement end of the interferometer electronics 
was described under the experimental section in Chapter IV. 
The analog bipolar signal was adjusted to be within a t1i0 V 


range and was connected to the analog input of the analog 
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to digital convexs@ongecircuitry. 

start Pulse and Clock Control Circuitry: A schematic 
diagram of the start pulse and clock control circuitry is 
shown in Figure 76. The start pulse and the digital clock 
were all controlled @y computer software. First, the 
comparator that converted the analog white light into pulses 
was set at a threshold so that normally three pulses were 
generated. Three instead of one were utilized because it 
was very difficult to maintain the threshold level to give 
only one pulse. This train of pulses was then converted into 
a Single pulse by a monostable whose delay time was set to 
about ten milliseconds. This ensured that only one clock pulse 
would trigger the JK flip-flop which then started the data 
acquisition and processing. Refer to the timing diagram 
iniFPisure 75 for the exact isequence of events.' The data 
acquisition and processing would start if and only if the 
computer issued a GO command (VIA PB4) and the white light 
pulse occurred during the forward movement of the mirror 
as distinguished by the synchronous signal. At this time, 
the laser digital clock appeared at the output of the 
NAND gates and thus triggered the conversion cycle of the 
analog-to-digital converter. If the computer did not issue 
a GO command, a logic level of '0O' appeared at the output 


of the NAND B gate and the system would be at the waiting 
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stage. Whenever a signal was acquired and processed, the 
microcomputer would issue a CLEAR command through the VIA 
Poo 1/0 line ho iclean bhe.JKiflip=flop.tosbeaseady for the 
next signal. 

Analog to Digital Conversion Circuitry: The analog to 
dieitali convenszon GireuLtry §is shownfin FPigure!77. It 
consists of a 12-bit compatible sample and hold (Analog 
Devices, Sample and Hold SHA 1A) and a 10-bit analog-to- 
digital converter (Analog Devices, Analog-to-Digital 
Converter ADC 10Z-002). The settling time for the sample 
and hold is five microseconds and the conversion time for 


the analog-to-digital converter is twenty microseconds. 


They were wired according to the specifications recommended 


by Analog Devices (69,70). 

The analog-to-digital converter operated under its 
internal clock during conversion of bipolar signals within 
+10 (V range.. It was not calibrated against a standard 
reference voltage nor zero adjusted for this preliminary 
investigation. Obviously, a 12-bit analog-to-digital 
converter which has been calibrated and zero adjusted is 
a better choice for the system, but none were available 
in the laboratory for this applicatiory. 

When the laser digital clock was used as the convert 


command for the analog-to-digital converter, the microcom- 
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pute KkepWtrachvot thesendeof.conversion.. This. was 
aecompiashed by the connection of the VIA CA1 line to the 
STATUS pin available in the analog-to-digital converter. 
In this way, the computer would generate appropriate 
correlation masks and control signals for data processing 
acrineend rofveviery feonversion. The .tO-bit .digital output 
in two's complement format, corresponding to the analog 
Signal, was continuously available to the data register X 
of the TRW multiplier-accumulator after every conversion. 
TRW Multiplier-Accumulator: The 64-pin multiplier- 
accumula terachip, MAC (LRW, Inc, TDC: 1009J.) is a high 
speed TTL LSI device capable of performing a multiplication 
between two parallel 12-bit data and product accumulation ) 
tntOsS Nanosccondsee7!) i.wlieeonsists.of two.12-bit data 
brepnustens,;pteand Yeebata, are doadedsinto.the registers x 
and Yiatether risang edge»of the. clockspulses, applied at 
CLK X»and' CLK Y respectively. Multiplication is performed 
between the two numbers in the data registers, followed 
byeproductsaccumulation at the rising edge of CLK P. The 
27= bat Outpuiaresister as tristate. controlled. It can be 
outputted in three segments, least significant 12-bit 


Co most significant 12-bit (Py 5-P5.) and extended 


14)? 
most significant 3-bit Grin P oe They are controlled by 


tie, is SMeandulsttristate controlesignals respectively. 
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A logic—-0*-at-2SL, TSM;-TSX-will-enabre the output 
drivers. 

ihe detailed, contigurataon of the mMAC is illustrated 
in Figure 78. The numerical system was configured in two's 
complement to be compatible with the analog-to-digital 
converter output. Since only 10-bit data were available 
from the analog-to-digital converter, the least 
significant two bits (X) and X,) of the data register X 
were grounded. When the end of conversion signal was 
detected by the computer, a correlation mask data point 
was fed into the data register Y of the MAC through 2 
VIA I10 lines (PBO and PB1). The reason why only 2 lines 
were used for a 12-bit correlation mask will be explained 
later along with the computer software. Once every- 
thing wast ready,” the *pers pheral’ controltline, CB2 or 
the microcomputer would generate the data loading clock 
pulses (CLK X and CLK Y) so that data were entered into 
the registers. The peripheral control line CA2 would 
then generate the multiply command (CLK P) to finish off 
the multiplication and accumulation process for one 
data point. The whole cycle was mepeaved for another data 
POL Until finished, 

Before data acquisition, it was necessary to 


preload the accumulator of the MAC with zeros. This was 
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achieved by the use of the non-accumulating mode of the 
MeO lozie “0* at/ACC). Data rezister Y was filled 

with zeros through the use of PBO and PBl I/0 lines. 
After the generation of the clock pulses (CLK X,Y, and 
CLK P), a zero product resulted in the-accumulator. At 
this point, the MAC was configured back to its accumulating 
Mode (logic ‘2' at ACC)’ to be. readyforiidatamacaui sition. 
The 27-bit accumulator output drivers were opened at the 
end of One scan by thesusepoft TSX, TSM, ZSL. command 
Signals. The accumulator content was loaded into the 
Computer through thegmultirplexing circuitry.) This was 
necessary because the computer ‘could only handle 

8-bit data at one time. 

Data Multiplexing Circuitry:s.«.he. MAC. output resister 
(27-bit) was multiplexed into the computer as three 8-bit 
and one 3-bit data points. A block diagram of the data 
me bor ekme: CAUCCUnbrye Vercnowl 1h) Preurem9 .) Hib 
consists of a 2-line to 4-line decoder and an array of 
Bf trisvaves. The controls of the tristates were soverned 
by the four output lines of the 2-line to 4-line decoder. 
The inputs to the 2-line to 4-line decoder were under 
computer control. A maximum of eight tristates were 
enabled at one time. The detail configuration of the 


2-line to 4-line decoder which was made up of six NAND 


nips, 
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Multiplier/Accumulator Output Control 


Control One 


Control Two 


2-Line to 
4-Line 
Decoder 


Rockwell 
AIM 65 
Microcomputer 


Controls 


Tristates 


8-Bit Data 


Multiplier/Accumulator Output Register 


FIGURE 79. Block diagram of the data multiplexing circuitry. 
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gates is shown in Figure 80. The PB6 and PB7 1/0 

lines were the two control lines from the computer. Only 
Ceenot tnese Tourjoutputs of the Z2-line to 4-line decoder 
would be at a logic '0' which would open up to a maximum 
of eight.tristates-at one time. The data transfer from 
the; MAC to the Port A of the computer is better understood 
with théewaid of Figure 81. During the transfer of data, 
the output data register of the MAC was enabled and this 
was easily accomplished by feeding a logic ‘'0" to the 
wiiepates Controls (Tob, Teja TSxX) of the MAC through 

the |FB5 line of the VIA, 

Computer Software: The program was written in BASIC 
and ASSEMBLER languages. It was initiated under BASIC 
and would request a wavelength input for the calculation 
of correlation masks. Correlation masks were calculated 
in vre=merm of square syntheuLc interferograms as 
described in Chapter IV and thus consisted of values ¢1, 
OFand sivonly. “Wherlc-bi® binary “equivalences or. these 
numbers are shown in Table 9. As it can be seen, 
the information of the correlation masks can be transferred 
with the aid of only two 1/0 lines. The PBO line of the 
VIA was used for the least significant bit and the PBl 
line of the VIA was used for the most. significant 11 bits 


of the coprelation masks to be transferred from the computer 
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Rockwell AIM 65 Microcomputer 


R6522 VIA Port B 
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FIGURE 80. Two-line to four-line decoder of the data 


multiplexing circuitry. 
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Rockwell AIM 65 Microcomputer 
R6522 Versatile interface Adapter 


PAO PAi PA2 PA3 PAS PAS PAB PAT 
[e) Oo Oo 


‘ort B I/O Lines 
PBS PRE 
© 0 OO 6 0 O 6 


Muttiplier/ 
Accumulator 
Output 
Contro! 


Pi70 
P220 Ai 
ee Baal’ 


TRW Multiplier/Accumulator Output Register (27-Bit) 


FIGURE 81. Multiplexing circuitry between the AIM 65 micro- 


computer and the TRW multiplier-accumulator. 
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to the MAC. The BASIC program stored the correlation 

masks (both sine and cosine) in decimal as shown in Table 
10. The decimal equivalence also included the information 
for pins PB2 to PB7. These were the control signals used 
to operate the system. The memory location of the 
correlation masks are indicated in the memory map shown 

ioe ave Lu, 

The BASIC program then requested the ASSEMBLER 
Subroutine which was entered into the computer with the 
aid of the Text Editor. The ASSEMBLER subroutine 
controlled the data acquistion and correlation data 
processing after preloading the accumulator of the MAC 
with zeros. The whole procedure was repeated twice, one 
for the cosine correlation and one for the sine correlation. 
The accumulator outputs stored in the memory by the 
ASSEMBLER subroutine were converted into the decimal 
mode by the BASIC program and outputed onto the display 
and printers’ Thetwhole precess could be repeated for 
another set of data with the same or a new correlation 
mask. This program was initially set for acquiring 
256-point interferograms and signal averaging capability 
was not available. However, the program can be modified 
to expand its capability. 


Program flow charts and listings are included in 
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TABLE 11. Memory map of the AIM 65 microcomputer. 


Type | Hexidecimal Decimal 


Bere aes Re aa a a a eee 


BASIC (internal use) 
BASIC program 
Assembler program 
Text editor 
Correlation mask (cosine) 
Correlation mask (sine) 
Results 
Index (used by Assembler) 


Assembler symbol table 
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Appendix III. For further details, the AIM 65 manuals 


(72,73) should be consulted. 


C. Evaluation of the Rockwell AIM 65 / TRW MAC Data 
Processing System 

Qualitative Aspects: The qualitative aspect of the 
system was tested out using a relatively simple signal, 
a He-Ne laser. The interferometric response was limited 
to the near IR-visible region because a silicon diode 
detector was used. The interferogram, which was a very 
Simple cosine wave carrying the characteristic information 
about the laser was correlated with masks generated from 
the wavelengths indicated in Table 12. The lines for 
the alkali metals were used because they are the prominent 
lines present in this region. The resulting correlation 
data based on two repetitions of the signal are shown 
in Figure 82. It can be seen that the correlation mask 
corresponding to the wavelength of the He-Ne laser 
shows a relatively high correlation result while all the 
others cornet. 

The second test involved the use of the signal 
from a magnesium hollow cathode lamp. The 
spectrum of the magnesium hollow cathode lamp is shown in 


Figure 83. This spectrum was obtained after Fourier 
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TABLE 12. Some spectral lines in the near IR-visible region. 


Spectral Line 


Code Number Wavelength (nm) 
i ee a en ee ee Pel RE ce Tee 


Na 589.0 
2 589.6 
ae 


1 766.5 
2 769.9 
1 780.0 
2 794.8 
x 652.1 
2 894.4 
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transformation of a 4096-point interferogram with the 

aid of the minicomputer used in the system described in 
Chapter IV. One has to keep in mind that the resolution 
of the signal obtained using the AIM 65 microcomputer 
(256-point interferogram) is not as good as the one 
obtained using the minicomputer (4096-point interferogram). 
The spectrum is shown here to indicate what spectral 
information is available from the magnesium 

hollow cathode lamp. From the spectrum, it can 

be seen that Mg 285.2 nm line is the prominant 

line present in the signal, followed by Mg 279.6 nm line 
and 280.3 line. This signal was obtained using the 

solar blind optical filter and the visible PMT detector. 
Correlation masks based on the wavelengths shown in Table 
6 were generated and used to perform real time data 
acquistion / processing with this system. The correlation 
data are shown in Figure 84. There is no distinct 
indication that Mg is the major component because the 
wavelengths indicated in Table 6 correspond to the major 
emission lines available from the ICP not the hollow 
cathode lamp. With the signal obtained from the hollow 
cathode lamp,the 285.2 nm line instead of the 279.6 nm 
line is the prominent line. This, together with the lower 
resolution offered by the system accounts for the high 
"background" for the other elements. However, the 


secondary correlation process was able to sort out 
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the information present in the signal. The 

elements that showed relatively high correlation results, 
namely dysprosium, lanthanum, magnesium, scandium, tin 
and thallium were further correlated at secondary ICP 
lines. Dysprosium does not have any other intense lines 
present in this region and could not be further correlated. 
fre the case of Banthanum;, it is very unlikely #that it was 
present due to the poor response of the detector for 

any spectral information outside the 240 nm - 360 nm 
region. The remaining six elements’ were then 

further tested with the secondary correlation analysis. 
The wavelengths used to generate the secondary cross- 
correlation masks can be found in Table 13. The secondary 
correlation results, as)shown in Figure 85, clearly 
indicate the presence of magnesium information in the 
Signal. From this example, the importance of the right 
choice of the spectral wavelengths used to generate cross- 
correlation masks is clear. 

Quantitative Aspects: Linear working curves for He-Ne 
laser (632.8 nm) and magnesium (285.2 nm) signals are illustrated 
in Figures 86 and 87 respectively. These “analytical 
curves" were obtained using neutral density filters to 
alter the signal levels. Thus, the concentration axis 


in both cases are relative in scale and have no resemblence 
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TABLE 13. Most sensitive ICP spectral lines of magnesium, 
scandium, tin and thallium within the 240-360 nm 


window. 


Spectral Line 


Element 
Code Number Beisaiieah?* hae (nm) 
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480 


He-Ne Laser (632.8 nm) 


420 
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Correlation Data (X103) 


120 


60 


fe) 0.2 0.4 0.6 0.8 1.0 
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FIGURE 86. “Analytical working curve" for the He-Ne laser 


(632.8 nm). 
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Magnesium (285.2 nm) 
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FIGURE 87. "Analytical working curve" for magnesium (285.2 nm). 
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at all to the actual level of the signals. These pseudo 
analytical curves are included to demonstate the 
quantitative aspect of this real-time data acquistion / 
processing system. 

The precision of the data obtained with this system 
was typically in the range of 9 to 10 % relative standard 
deviation (rsd). It can at best attain a rsd of 4 to 5% 
depending on the wavelength and the amplitude of the 
Signal. This modest level of precision is limited by 
the fact that phase correction has not been performed on 
the same set of data. Any change in phase or amplitude 
between scans would contribute to the scattering of 
the readings. The solution to this problem is to implement 
phase correction on a single scan of signal with cosine 
and sine correlations performed simultaneously. This 


can be easily achieved with the use of two MAC's. 


DasCeonclusionsvand gProspects 

At the present stage, the Rockwell AIM 65 micro- 
computer / TRW MAC system cannot function properly during 
thepoperation lof ithe eheP, duectorradboirequency minteri erence 
originating in the plasma power supply. This is the 
reason why no actual ICP experimental data could be 


presented. This interference problem should be able to 
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be completely eliminated by properly packaging the whole 
system and also by the use of by-pass capacitors on every 
Single electronic chip that has been used in this system. 
Further details on the use of by-pass capacitors can 
be obtained in reference 74. These have not been tried 
at this stage because the system requires a major 
modification to incorporate two MAC's into the system for 
better quantitative data. 

All the results obtained are all based on a 256-point 
interferogram and the system can certainly be expanded 
to increase the resolution (i.e. decrease spectral 
interferences). There is actually no limitation on the 
length of interferograms that the system can handle, 
because the signal obtained is processed in real time 
and is not stored in the memory of the computer which 
normally imposes a limitation on the number of data points 
to be processed. However, in the case of long interferograms, 
the accumulator register of the MAC has to be checked 
for overi low. 

The data processing system can be modified to a 
single channel slew scan type processor with the use of 
2 MAC's or a multichannel direct reading type processor 
with the use of more than 2 MAC's. In the case of the 


latter design, one pair of MAC's would be dedicated for one 
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spectral line and true simultaneous multielement analysis 
can be accomplished. The use of the square synthetic 
interferogram correlation masks can be very beneficial to 
the system. The correlation masks have values of +1 and 
-l with rare occasions that the value is exactly equal 
to zero. Thus, the zeros can be replaced by either +1 or 
-l1 value with no substantial effect on the selectivity of 
the cross-correlation masks. A significant saving of 
computer memory for storing correlation masks can be 
achieved. High speed real time correlation analysis is 
possible if all the masks can be stored in ROM, instead 
of wasting computation time for the generation of the 
correlation masks. Moreover, correlation masks can be 
input into the MAC's more efficiently because one 
byte (8-bit) of data output will be able to feed 8 MAC's 
with correlation masks simultaneously. The AIM 65 
microcomputer is limited by its 4K RAM and it should 
be expanded to increase its versatility. A large number 
of correlation masks can then be stored in ROM and recalled 
whenever they are necessary. 

Based on these, a versatile, flexible and inexpensive 
real time correlation based data processor should be able 
to emerge for the I1CP-Fourier transform spectrometer 


system. In fact, it should be treated as a general 
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purpose data processor for interferometric signals. 
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CHAPTER VI 
Summary 


In this study, the feasibility of automatically 
processing spectral data by cross-correlation techniques 
has. been investigated. In particular, cross-—correlation 
functions at t=O point have been utilized. The basic 
approach involves cross-correlation of the raw spectral 
Signal with a noise-free mask of the sought-for spectral 
pattern at zero phase shift. Various types of cross- 
correlation masks have been evaluated as to their 
effectiveness of spectral signal detection. These cross- 
correlation masks were either experimentally or computer 
software generated. 

For an inductively coupled plasma - photodiode 
array spectrometer (ICP-PDA) system, several analog 
and binary spectral masks originating from high signal-to- 
noise ratio spectra of the sought-for elements were 
utilized to process the analytical information available. 
Both analog and binary cross-correlation masks clearly 
showed their effectiveness in extracting sought-for spectral 
information from complex multielement emission spectra. 
Spectral overlap problems were overcome by the use of 
highly selective cross-correlation masks. Furthermore, 
with the use of these cross-correlation masks, all the 


spectral information related to the sought-for chemical 
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species available from the system were fully utilized. 
In this way, the analytical data obtained were better than 
those obtained using the conventional single peak methods 
in terms of sensitivity and precision. Due to the limited 
spectral coverage of the ICP-PDA system, correlation 
analysis could not be utilized to its full potential, but 
the prospects of this technique certainly look good if 
photodiode array spectrometers, or in general, spectrometers 
based on image detectors ‘can be developed with wide spectral 
coverage. 

In the case of an inductively coupled plasma - 
Fourier transform spectrometer (ICP-FT) system, due to its 
inherent property of precise wavenumber calibration, 
analog and- square cross-correlation masks were computer 
generated for use in the correlation analysis of interferograms. 
This, together with the wide spectral coverage provided by the 
Fourier transform spectrometer provide a high degree of 
flexibility in the use ‘of the system, The cross-correlation 
masks are highly effective in obtaining quantitative and 
qualitative information directly from the interferograms 
obtained by the spectrometer. In carrying out elemental 
identification on the signals, a preliminary correlation 
operation is performed for all the elements. This procedure 


essentially shortens the list of elements that are suspected 
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to be present in the sample. Further investigations by 
correlation analysis on the suspected elements yield the 
final analytical results. A consequence of the simplicity 
of this correlation method, as compared to the usual FFT 
data processing method, is that a real time data processing 
system using a low cost microcomputer could be developed. 
The system utilizes a hardware multiplier and accumulator 
to carry out the actual correlation computation. The 
preliminary results obtained look promising and this 
approach warrants further investigations. 

The comparable capability of the simplified 
cross-correlation masks, namely the binary masks for the 
ICP-PDA system and the square masks for the ICP-FT 
system when compared to the analog cross-correlation 
masks has certain advantages when implementing the 
correlation techniques. The correlation procedure 
becomes Simple and fast, and binary and square masks 
are simple to store in read-only memories (ROM). 

This is important because storage of cross-correlation masks 
is necessary for high speed data processing. 

In conclusion, spectral pattern recognition tasks 
are among the most difficult and complex operations to 
program and implement with electronic and computing 


technology. In fact we have some way to go before we 
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can rival the combination of the human brain and eye in 
such tasks. However, cross-correlation based procedures 

are beginning to provide a modest level of capability in 
automatic interpretation of complex data arrays and one 

can certainly foresee the important role of correlation 
techniques for data processing in simultaneous multielement 


analysis, or in general spectrochemical analysis. 
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FOOTNOTES 


The notation used to label the axes throughout the thesis 
follows the version adopted in the computer graphics 
industry. For, example, in: Figure 3, the amplitude axis 


represents values ranging from 0 to 6000. 


In this particular example, the dynamic range problem 

is imposed by both the emission of the plasma background 
and the appearance of the noise spike. In cases where 
plasma background emission is low, such as in the solar 
blind spectral region, the dynamic range problem may be 


solely governed by the amplitude of the noise spike. 


The sodium lines appear between the lithium and potassium 


lines because of alaising (28). 
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APPENDIX I 


Alternate Methods of 


Recycling Fast Fourier Transform Algorithms 


In this appendix, several recycling implementations 
other than the one presented in Chapter II when Scheme 2 
is used as the primary FFT method will be illustrated. The 
computer software is identical to the one described 
in Chapter II and will not be repeated here. 

let us recallethat scneme 2 does produce zeal and 
imaginary output arrays and the imaginary array was ignored 
in the recycling method presented as Scheme 4 in Chapter II. 


It would seem that a recycling method as illustrated below: 


bars ft 
I (even points) eos iM 


1/25R 


(5) 


it -I 


should also be valid. In this scheme the real output array 


is put into the real input array of the second FFT and the 


negative of the imaginary output array is put into the 
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imaginary input array of the second FFT and the FFT 
implemented in the normal way. The real, imaginary and 
amplitude output arrays for the second FFT are shown in 
Figure 88. This method works, but now the data occupy 
only hailf of the output arrays. 

Another approach is to ignore the imaginary output of 
the first FFT and recycle using Scheme 1 to set up the 


second FFT. This is outlined below: 


aa (odd points) PFT R 
DATA ESD 
oe (even points) WITH SORT al 


. get 


ae TC ZEROS) 


(6) 


This also "works" as shown in Figure 89, but now the output 
arrays contain redundant information in the upper and lower 
halves. 

Finally, a pseudo double sided interferogram similar 
to the interferograms generated by Fourier transform 


spectrometers can be constructed from the real output of 
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FIGURE 88. (a) Real, (b) imaginary and (c) amplitude 
output arrays for Scheme 5 inverse FFT 
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FIGURE 89. (a) Real, (b) imaginary and (c) amplitude output 


arrays for Scheme 6 inverse FFT implementation. 
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the Scheme 2 FFT. It is simply necessary to "reflect" it 
about the origin generating a new array, twice the number 
of the original data points as shown in Figure 90. The 
real, imaginary and amplitude output arrays for the second 


FFT when processed as outlined below: 


R (odd points) R 
FFT 
DATA 
WITH SORT I 
I (even points) 
Reflect to 
pseudo 
interferogram 
R PPT R=sGodd= points) 
A C2) 
i WITH SORT il (even points) 


are shown in Figure 91. The amplitude array now contains the 
original data. 

While Schemes 5-7 do represent somewhat unconventional 
inverse FFT routes they do clarify the information 
content of the various input and output arrays and illustrate 
that no one method is solely acceptable. In addition Scheme 
7 illustrates the relationship of "double sided" interfero- 
gram type signals as generated by Fourier transform spectro- 
meters to the output arrays when conventional spectra are 


transformed to the Fourier domain. 
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FIGURE 91. (a) Real, (b). imaginary and (c) amplitude output 


arrays for Scheme 7 inverse FFT implementation. 
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several analogous manipulations of the real and 
imaginary output arrays for recycling when Scheme 1 was 
used to implement the first transform were investigated. 
None proved useful. In general, they generated severely 


overlapped information. 
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APPENDIX II 
Modifications of the Michelson Interferometer 


The Michelson interferometers in our research labor- 
atory have evolved over more than 15 years. The original 
version of the interferometer was designed and built by 
Chaney (29,30) at the University of Michigan and was later 
adapted by Horlick (31) in 1969 at the University of Illinois. 
This design was then significantly modified into the 
present form by Yuen (32) and Hall (33). In this appendix, 
modifications that have been made to the interferometers 


Since then, up to the time of writing will be described. 


A. Photomultiplier Tube Dynode Chain Circuitry 

Two changes have been made on the dynode chain 
elecurttry. Bie l05emvenerarad capacitor was originally 
wired to the negative voltage supply (-HV) by mistake. 
The configuration has now been changed so 
that the capacitor is connected to the ground. This will 
ensure all the high frequency noise that is going through 
the dynode chain to be grounded. The other change is the 
conversion of the voltage output circuitry to the present 
current output configuration. This simply involves the 
removal of the load resistor and enables the current 
amplifier to be described later in this appendix to be 


utilized. 
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The final configuration of the photomultiplier tube 


dynode chain circuitry is shown schematically in Figure 92. 


B. Photomultiplier Tube Housing 
A 0.5 inch diameter UV grade quartz lens (Oriel 

Ai1-621-08) of focal length 13 mm is inserted between the 
entrance aperture of the PMT housing and the PMT. It is 
placed at’a distance equal to its focal length from the 
aperture so that there is a significant increase in the 
amount of light that will reach the PMI. No experimental 
test has been carried out to see what is the difference 


between signal levels obtained with and without the lens. 


C.. Stlicon Optical Detector 


A new sidicon dvodewdetector from 2R Industries, Inc. 


(catalog number: 8016L) is now being utilized. This detector 


provides better sensitivity than the old silicon diode 
because it has a large sensitive area and is equipped 
Withéa lens fomrfoctisamp@wlight onto the silicon dicde, 
Again, no test has been done to check on the performance 


of these detectors. 


D. Optical Filter Mount for the Detectors 


Optical filters are invaluable for resolving dynamic 
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range problem imposed by the multiplex property of the 
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Michelson interferometer (see Chapter IV). A mount designed 


for optical filters is attached to the detector mount. 


Optical filters such as those obtained from Oriel 


Corporation of America can be easily put in place in front 


of the detector. The. optical wmount,.can take anya inch 


diameter filters with thickness varying up to approximately 


tan). 


E. Micrometers on the Fixed Mirror. Mount 

The Michelson interferometer consists of a fixed 
mirror and a movable one. Final signal optimization is 
normally. achieved by the alignment of, the. fixed mirror. 
In the case of UV signals, the alignment of this mirror 
is crucial and is not easy to achieve using a 
screw driver. Two differential micrometers obtained from 
Lansing Research Corporation, Ithaca, New York (model 
number: 22 505) are incorporated into the fixed mirror 
mount. Perfect alignment of this mirror can now be 


achieved very easily with the use of these micrometers. 


F. Photomultiplier Tube Power Supply 
The PMT is now powered by a high voltage supply 


obtained from Keithley Instruments (model number: 244). 
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The voltage can vary from -200 V to -2200 V. 


G. Signal Measurement Electronics 

The signal measurement channel has been redesigned with 
new electronic equipment . The current output from the 
detector is amplified with a Keithley 427 current amplifier 
before it is electronically filtered with a Krohn-Hite 
3343 filter. The current amplifier can have a gain from 


107 to 1021 


with d.c. suppression. Additional amplification 
can be obtained from the electronic filter (maximum 40 GB). 
This is usually configured as a bandpass filter. Final 
adjustment on the gain of the signal is provided by the 


PAR 225 low noise amplifier before the signal is input 


to the computer. 


H. Old Mirror Drive System 

The electrical wires that provide the power to move 
the mirror have been changed to smaller wires to provide 
smoother drive. In addition, some stop pins have been 


inserted into the mount to secure the springs in place. 


I. New Mirror Drive System 
Perhaps, the major modification to the interferometer 


system is at the mirror drive system. Two interferometer 
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systems are currently available in the laboratory. One 
with the old mirror drive system and the other with a 
new mirror drive system that is going to be described in 
Rais seculon. 

A linear actuator (or linear moving coil motor) 
purchased from Kimeo, Inc., San Marcos, California is 
used to drive the mirror in one of the two interferometer 
systems. This drive is capable of delivering a peak force 
of 8 pounds. With this new drive system, preliminary data 
show excellent stability even without the electronic 
phase locked loop and is certainly a great improvement over 


tie. OolG mirror drive. 
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APPENDIX III 
Computer Software 


he Program Flow (Chart agd Listing for the Fast Fourier 
Transformation Program Based on the Cooley-Tukey 


Algorithm 
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START 
INPUT NUMBER 
OF POINTS 
INPUT DATA 

FROM DISK 


RECONSTRUCTION 


OUTPUT REAL 
AND IMAGINARY 
ARRAYS ONTO 

DISK 


CALCULATION OF 
AMPLITUDE 
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OUTPUT 
AMPLITUDE 
SPECTRUM ONTO 
DISK 


FIGURE 93. Flow chart for the FORTRAN FFT program. 
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THIS PROGRAM WILL PERFORM THE FAST FOURIER 
TRANSFORMATION USING THE CQOLEY-TUKEY 
ALGORITHM 


FROGRAM IS EXECUTED UNTER A RT-11 OPERATING 
SYSTEM WITH A DEC FOF 11710 MINICOMFUTER 


SUBROUTINES REQUIFEL? 
(1) FORLIE 
(2) SY Silene 


TIMENSION X£4996) 2Y64096) 21. (20) 


COMPUTER REQUESTS AN INFUT TO INDICATE WHETHEF 
THE FFT IS WITH QF WITHOUT THE SORTING ROUTINE 


eRe lis FE TAROUT INE 
INFUT ONE DATA ARRAY 
- OF0 POINTS INTO REAL ARRAY X CN) 
- EVEN FOINTS INTQ IMAGINARY ARRAY YON) 
“MERESPRQEALEER WITH ZEROS 
per ae con 
—~w RECONSTRUCTION 
- TATA OUTFUT 


NON-SORTING FFT ROUTINE 
- INFUT TWO DATA ARRAYS 
~ FFT 

- DATA OUTFUT 


WRITE 72100) 

FORMATC’ “%s’O=SQRT 1=NO SORT“? 
READ CS 2120). LSORT 

FORMAT (I2) 


COMPUTERMREQUESTS INFUT OF NUMBER OF FFT POINTS 
IN THE FORM OF 2¥*N 


MAXIMUM NUMBER OF DATA FOINTS IS 49094 (N=12) 
FOR THIS FROGRAM 

LerMaTAgeGONeseessUE TO THE SIZE OF THE “COMPUTER 
MEMORY 


WRITE C7140) 

FORMAT(’ “’s“INFUT N (2*#*N) 7%) 
READCS 2120) N2FOW 

NF =2**N2FOW 


DATA INFUT 


ONE INFUT ARRAY FOR SORTING ROUTINE 
TWO INFUT ARRAYS FOR NON-SORTING ROUTINE 
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CALL ASSIGN(192 “DKS DATA.LAT’ s-1s OLD’ ’NC/% 9d) 
READICLO) CX (IT) se T=1)NPF) 

IFC ISORT.EQ.0) GO TO 9 

REANICLOCYCI) es IT=1°NF) 

ENDE TLES LO 

CALL CLOSE (10) 


SORTING OF DATA ARRAY FOR SORTING FFT ROUTINE 


Ort DATA FOINTS INTO REAL ARRAY 
EVEN DATA FOINTS INTO IMAGINARY ARRAY 
BESTIEILEEM WITH ZEROS 


LRGLoUn ets 2 GO TO 12 
NF T=NF 7/2 

N=NF* 

UG 2 Deal 

X(T) =XCIkK2-1) 
YCIV=XCI¥2) 

NO 3 I=NFT+1 NF 
Y(I)=0,.0 
X(T)=0,0 


OOOO moO om oom mF 
COOLEY-TUKEY FFT ALGORITHM 
OOOO OOOO OOK 


NTHFOW=2¥¥NOQFOW 
N4FOU=N2FOWZ2 

IF (N4FOW) 69263760 

NO 61 IFASS=1-:N4F OW 
NXTLTH=2*¥ (N2FOW-2¥ITF ASS) 
LENGTH=4¥NXTLTH 

SCALE=6. 28416S53/FL0ATCLENSTH) 
rd 61 J=1-°NXTLTH 
ARG=FLOAT (J-LI*SCALE 
Ci=COS¢CARG) 

S1=SINCARG) 
Ce=CieGl—-Si*tsi1 
S$2=C1i¥*S1+C1x*S1 
C3=CURC2L-Si¥S2 
S3=C2k514524C1 

ro 61 ISQLOC=LENCGTH:NTHPOWsLENGTH 
Ji=ISQLOQC-LENGTHt+J 
J2=JL+NXTLTH 
J3=J2+NXTLTH 
J4=J3+NXTLTH 
Ri=X¢J194+X¢C IS) 
R2=X%CJ1)-X¢CI9) 
R3=X(J2)+X¢(JS4) 
R4=X%€J2)-X¢(J4) 
FIL=Y¥¢( JL) tY (JSS) 

Pi2=¥ (3021S) 
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FIS=Y(J2)4Y(14) 
FI4=Y¥(J2)-YC54) 

XC JLI=RLI+tRS 

¥ O41 Y=Rl bP IS 

IF (J-1)64962964 

MC ISZVECLIECR2¢F 14491 ¥ CF T2-F4) 
YC 435) =-SL¥(CR24+F 14) 4C1¥ (CF I12-k4) 
X% CA? eC2¥ CRA RSRS SS OR TIF TS) 
Y GID =-S2eGnIoRS+E2xk (Fr Li-F IS) 
XC 14) =C3¥KCR2-FI4 +S 3k (R44FI2) 
Y(J4)=-S3¥ (R2-FI4)+C3k(R44F1I2) 
COON OF 

XC IZSI=RLI+F 14 

YC IZ)=FI2-R4 

¥(J2)=R1-R3 

Y GID Shi -Pie 

X(1J4)=R2-FI4 

YC4I4)=R4+FI2 

CONTINUE 

IF CN2SPOW-2¥N4FOW) 4652462465 

TQ 67 J=L»NTHFOWs2 
RI=XCJtXCIS+1) 

Rook C= 7st 1) 

FI HVC») 

Fe Sy) e) 


Y QISFId 

XC I¢F1)=R2 

Y GEEPY=PI2 

Na 68 J=1913 

L¢so=1 
IPCI-N2FOMOGSt6o » 68 
L¢ I) =2e# CNQFOWFI1- J) 
CONTINUE 

T=] 

Lr=C¢t3) 

NTHFOW=L 012) 
ISQLOC=L(11) 

J=L¢10) 

N2OFL=L (9) 

N2=L(8) 

NF 2MJ=b 07) 

L8=L. 06) 

N2FOW=L (05) 

NAR OWL €4 ) 
LENGTH=L¢3) 
NXTLTH=L¢2) 
IFASS=L(1) 

Vas Gor Itty 

NO 6901 J2=JiysNTHFOWsL 1 
0 601 J3=J2*»ISQLOC:NTHFOU 
nO 601 J4=J3»J2TS@LQC 
N10 601 JS=J4eN2F 1+ LI 
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Jé=I5 9N29N2P4 
J17=159NP2MI-N2 
J8=J7+L8+NF2MJ 
JI9=JI8»NOFOWrLE 
11019 NAP OU »N2FOU 
Jil=J10:LENGTHsN4FQu 
Ji2=J11.NXTLTH: LENGTH 
JI=J12:IFASS:NXTLTH 


troC La Tele els 601 


R=XCIJ) 


bia Nepal eo ee ALA) 


XC IT) =F 


Pay Cl 
VCL I= Y Caio 
ACS ING Ea) a 


IJ=ITJI+1 


SOO OSS i kk ip 
COOLEY-TUKEY FFT ALGORITHM ENDS 
SOO oi ob om Om fofabak tok 


RECONSTRUCTION FOR THE SOFRTI 


a) iad Ct UES Sg ae co er) Mel C1 Bp M8 Wg 
4 


AFIS ss 


:1L415927/FLOAT ON) 


C1=COSCARG) 
S1=-SINCARG) 


CiJjx=1, 
S1IX=0, 

N2=N/2 

NOP IEN2+1 
Je2eN2r i 


nd 
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NF 2MJENt2~J 
SORFRL=X(1)4+X¢NF 2M) 
SORLE=Y C0] CNE SM) 


Rise 1 IX 


C1JIX=C1IZX¥C1-S1IX*S1 
SiJX=R¥51+S1IX¥C1 
SORR2Q=*¥ CJ) -X CNP 2M 1) 

SOR Leo) (oa) ON aii 
SOQRFR3=C1IX¥SQRR2-S1I*X*SORI2 
SORT 3=C1IX*SORT2+S1IX¥SORR2 
Y(J)=0.5* (SORTI-SORKS > 

XC J=0,5e (SORRI+SORI 3S) 
Lela 2 Eel Osyeieee Os7 2 

¥ CONF 2M) 0-48 -SORTI-SORE 3) 
X(NF2MJ)=0,.5¥* CSORRI-SORIS) 
CONTINUE 

X(1)=XCLI+FVC1) 


Y(1)=9, 


VATA OUTFUT ONTO THE DISK 


REAL ARRAY OUTFUT UNDER THE 


NG FFT FOUTINE 


FILE NAME 


244, 


eal’? 


7 7 a ) ; ‘ ; = ms 
" : 7 ‘ 
113 ann | ‘on - 
a iy ena — 


cSMnet, at toa. Ay cel 
S4Veabs=tk - toa 
UPC $8.46 TL eOl 1 PO8 7 
Gieateasce sat £0e og 

WORM sSHOTAM.e) Ott : 50a Od: 
Fhe ere hey 208: “94 
PMP Lenk tOerad™ |, 
AION SLL He 08-00 | 
AOS1O L435 ett rlIP . 
LE)X=a 


7 
7 
' 


; (TL yee (k T3% | 

=e ‘et IL) 
(CL ¥<IF 

| (LEPPSCL ID | 
ct } (RISC 
We, 2 PRM Ra b Tia Oi 
. 7 
Hee RNa rKrs se reeenes #, 
24 i hey Rae Ta9 YANUT "3 Oe: . 
eRe ee rua Oey her eer? *¥ CREF EES ERE = mn 


sary — 
, —- 


S4T Tae ay aus MOL tua ena aay 
i OF OP 1.05, Tange /aeb. on 


SRA a a / 
aX 4 ro - . Y = 
(oa! NaI \4 p eee > 25 ‘aad aa on rs 

- 


tS ie 


j “ rete > 
; | a; | BN af PS. 
7 _. Seer vA . 


ibis tek os es 


t~ erhek a 
, CLM) X FOE) 0c 
| INS YO ee 
) 7 i, eee 
We bed AE Ciessieh onal 
. 7 Eh at ah SEE 
, he . CUMESA Se 1L Ve cee 


| os CLMSSAYY + | 4) 
ERROR SRE Be SAnce NL Z 738 
SRS eM S95 12080 
Ai « 2: e- EIAGE key hes 


oui 02h. Oa : 
7 ; a ee ove oat 
hee -FIANe =e, 3 a MEA) Y 

ME 20 at , a Per hea eY . 


rT N 
> a 


mo Leal (ea) Leah lee 


(eal (63) Le 


fC) C70 OCG Gt Gi) 


245 


—~REAL «DAT 
IMAGINARY ARRAY OUTFUT UNDER THE FILE NAME 
-IMAG. DAT 


CALL ASSIGN(10*“DKtREAL DAT’ s11+s’NEW’s“NC’% 1) 
WRITE CIO) (X¢€T)sT=1+NF) 

ENDETLE 20 

Unit, CLOSE CLG) 

CALL ASSIGN(10: “DN? IMAG. DAT’ sii: “NEM © /NC’%»1) 
WEITER LOD CY CLs belo Nr) 

ENDER ILE «10 

CALt. CLOSE Gis 


CALCULATION OF THE AMPLITUDE SPECTRUM 


POP Sti is he Ne 
XCTI=SQRTCXCTIEXCTIFYCIKYCT)) 
CONTINUE 


QNFPLITUNE SPECTRUM OCUTFUT ONTO> THE SUIS UNDE 
THE ORIEE NAVE 
-AMPL «DAT 


TALL ASSIGN(10+°UR TQMPL DAT! sills ’NEW os ONO’ ed) 
WRITE CLO) ¢X¢CI)2T=1l9NF) 

PNUP TCE LO 

CALlo -CLeSe cI) 


FROGRAM ENDS 
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ROGEE a ioe ING 
ROK OOK Ke 


SUFERVISOR? DR. GARY HORLICK 


JUNE 1780 

CHEMISTRY DEPARTMENT 
UNIVERSITY QF ALBERTA 
EVIMONTON» ALBERTA 
CANATIA 
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B. Computer Software for the AIM 65 / TRW Multiplier - 
Accumulator Data Processing System 
(a) Program flow chart and listing for the main program 


in BASIC. 


ft | 
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START 
INPUT 
WAVELENGTH 


CALCULATION AND 
STORAGE OF 
CORRELATION 


DATA ACQUISITION 
AND CORRELATION 
(ASSEMBLER 
SUBROUTINE) 


CONVERSION OF 
BINARY RESULTS 
INTO DECIMAL 


OUTPUT 
CORRELATION 
DATA 


NEW 
CORRELATION 


FIGURE 94. Flow chart for the BASIC program. 


10 
LS 
20 
25 
30 
40 
50 
60 
70 
80 
90 
100 
0 
120 
230 
140 
150 
TSS, 
160 
170 
7S 
180 
1o0 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 


INPUT "WAVELENGTH ,NM";W 
W=3975 .99/W 

FOR I=0 TO 255 
A=I*W 

R=COS (A) 

K=0 

IF R<O THEN J=255 
IF R=0 THEN J=252 
IF R>O THEN J=253 
IF K=l THEN 130 
POKE 3072+1,J 
Meek 

R=SIN (A) 

GOTO 50 

POKE 3328+1,J 
NEXT 

POKE 4,0:POKE 5,7 
PRINT "WAITING" 


GET AS 

IF AS="" THEN GOTO 160 
PRINDEDE PATIENT" 
DUM=USR(T) 

I=0 

IF I=0 THEN GOTO 230 
G2=Cl 

I=4 


A=PEEK (3584+T) 
B=PEEK (3585+I) 
C=PEEK (35 86+TI) 
D=PEEK (3587+TI) 


C1=D*1.6777216E07+C*6 .5536E04+B*2 .S6E02+A 
IF Cl>6.7108864E07 THEN C1l=C1-1.34217728E08 


IF I=0 THEN 210 
RES=SOR(Cl*E1+C2*C2) 
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2 ON 


te 


310 
320 
530 
340 
350 


PRINT! RES 

PRINT"O=NEW MASK 1=NEW DATA" 
TNPUDEE 

IF I=0 THEN 10 

GOTO 150 
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(bo) Program flow chart and listing for the subroutine 


program in ASSEMBLER. 


7? 


i 
ain 
snitvostiue- ont. tot tai bias Pato - weak satgout.. 
a) whe ssianeie aA “ar aaxqort 
va ay ib q 


ely 


MAIN PROGRAM 
IN BASIC 
PARAMETERS 
SET UP 


PRELOAD TRW 
ACCUMULATOR 
WITH ZERO'S 


CONVERSION 
FROM ADC 


OUTPUT 
CORRELATION MASK 
AND APPROPRIATE 
CONTROL SIGNALS 
AND CLOCK PULSES 


LOAD TRW OUTPUT 
REGISTER INTO 
AIM 65 MEMORIES 


RETURN TO 
BASICS 


FIRURE 95. Flow chart for the ASSEMBLER program. 


LOOP 


LONE 


*=$0700 


LDA 
STA 
LDX 
LDY 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 


#$00 

SOFO0 

#$00 

#SFF 

# SFF 

$A002 

#$00 

$A003 
#%11001100 
SA00C 
#%11100000 
$A000 
#%11101100 
SA00C 
#%11001100 
SAO0C 
#%11001110 
$A00C 
#%11001100 
SA00C 
#%11111100 
$A000 
SA0O0D 
#%00000010 
LONE 

$OFO0 

LTWO 
$0D00,X 
LTHR 
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; INDEX FOR SINE OR COSINE 

; MULTIPLICATION 

;X REGISTER AS POINTER 

;Y REGISTER AS NUMBER OF POINTS 
;PORT B FOR OUTPUT 


;PORT A FOR INPUT 
;CAl, CBl INPUT, FALLING EDGE MODE 
;CA2, CB2 OUTPUT, MANUEL PULSE MODE 


;PRELOAD MULTIPLIER WITH ZEROS 


sGENERATE CLOCK. X,Y ,.WLITH CB2 


; GENERATE CLOCK P WITH CA2 


; READY 
CHECK FOR END OF CONVERSION ON CA1 
;YES, GO AHEAD 


; LOAD CORRELATION MASK AND CONTROL 
;SIGNALS ONTO PORT B 
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| ‘Bouse Ave 


es 
Bi 


LTWO 
LTHR 


GPIV 
LSIX 


LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
INX 
DEY 
BNE 
LDA 
BEQ 
LDX 
JMP 
LDX 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 


$0C00,x 
$A000 
#%11101100 
SA00C 
#%11001100 
SA00C 
#%11001110 
SA00C 
#%11001100 
SA00C 
SA001 


LONE 

SOFO00 
LFIV 

#$04 

LSIX 

#$00 
#%11100100 
$A000 
#%11101100 
$A00C 
#%11001100 
$A00C 
#%11001100 
$A000 
$A001 
$0OE00,X 
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;GENERATE CLOCK X,Y AND P 


;CLEAR INTERRUPT FLAG 

; INCREMENT POINTER 

;DECREMENT COUNTER 

;REPEAT 2F NOT FINISHED 

;LOAD ACCUMULATOR OF MULTIPLIER 
;INTO AIM 65 


;LOAD MULTIPLIER WITH ZEROS TO 
; ENSURE NO FURTHER ACCUMULATION 


;LOAD PRODUCT IN 8-BIT SEGMENT 
;WITH LEAST SIGNIFICANT FIRST 


| ooLrekitey Aaa 
me i S08 ATA i 
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LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
AND 
STA 
LDA 
BNE 
INC 
JMP 
RTS 


#%01001100 
$A000 
SA001 
$0E01,X 
#%10001100 
$A000 
SA001 
$OE02,X 
#%00001100 
$A000 
SA001 
#%00000111 
$0E03,X 
SOF00 

LFOU 

SOFOO 

LOOP 


;CHECK IF COSINE MULTIPLICATION, 
;REPEAT WITH SINE MULTIPLICATION 


;RETURN TO BASIC 
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